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The differential cross sections for inclusive neutral pions as a function of transverse and longi¬ 
tudinal momentum in the very forward rapidity region have been measured at the Large Hadron 
Collider (LHC) with the Large Hadron Collider forward detector (LHCf) in proton-proton collisions 
at y/s = 2.76 and 7TeV and in proton-lead collisions at nucleon-nucleon center-of-mass energies of 
-ysNN = 5.02 TeV. Such differential cross sections in proton-proton collisions are compatible with the 
hypotheses of limiting fragmentation and Feynman scaling. Comparing proton-proton with proton- 
lead collisions, we find a sizable suppression of the production of neutral pions in the differential 
cross sections after subtraction of ultra-peripheral proton-lead collisions. This suppression corre- 
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sponds to the nuclear modification factor value of about 0.1-0.3. The experimental measurements 
presented in this paper provide a benchmark for the hadronic interaction Monte Carlo simulation 
codes that are used for the simulation of cosmic ray air showers. 

PACS numbers: 13.85.-t, 13.85.Tp, 24.85.+P 


I. INTRODUCTION 

Observations of high-energy cosmic rays with energy 
above 10^^ eV provide key information for a yet un¬ 
established origin(s) and acceleration mechanism(s) for 
cosmic rays. The compilation of current observations 
reveals kinks in the energy spectrum that agree with 
the turning points in the mass composition [1] at ^ 
3 X 10^^ eV (the so called ’knee’) and provide a consis¬ 
tent description of the transition from Galactic to ex- 
tragalactic cosmic rays at ^ 5 x 10^® eV (the so called 
’ankle’). In particular, a cut-off feature of ultrahigh- 
energy cosmic rays (UHECRs) at ~ 5 x 10^® eV is sup¬ 
posed to existence of Greisen-Zatsepin-Kuzmin HE] cut¬ 
off, while the source and propagation of the UHEGRs is 
still a mystery j?]. In order to grasp the experimental 
signature of the source of UHEGRs, and to understand 
consistent picture of transition from galactic component 
around 10^"^ eV, many extensive air-shower experiments, 
including on-going UHEGR observatories (i.e. Auger H 
and Telescope Array m) have collected the data on en¬ 
ergy spectrum, mass composition, and arrival direction 
of UHEGRs high energy cosmic rays over the past few 
decades |7]1|. 

It is important to note that critical parts of the 
analysis still depend on Monte Garlo (MG) simulations 
of air shower development that are sensitive to the 
choice of hadronic interaction models. Therefore differ¬ 
ent hadronic interaction models, which simultaneously 
predict the soft and hard QGD interactions, provide dif¬ 
ferent viewpoints even using exactly the same data com¬ 
pilation Hiini. Gurrently the lack of knowledge about 
forward particle production in hadronic collisions at high 
energy hinders the interpretation of observations of high- 
energy cosmic rays m- 

Here it should be remarked that the LHG at GERN has 
so far reached 13TeV centre-of-mass energy in proton- 
proton {p + p) collisions. This energy corresponds to 
the cosmic ray energy 9.0 x 10^® eV in the target rest 
frame which is well above the first turning point in the 
mass composition of primary cosmic rays from proton 
dominated to light nuclei dominated, namely the knee at 
approximately 3 x 10^® eV mi. The data provided by the 
LHG in the forward region, defined as the fragmentation 
region of a projectile particle, should thus provide a useful 
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benchmark for the MG simulation codes that are used for 
the simulation of air showers. 

The energy in the laboratory frame converted from 
the collision energy in p + p collisions at = 7 TeV 
(Aiab = 2.6 X 10^® eV) is two orders of magnitude lower 
than the ankle region where a transition from galactic to 
extragalactic cosmic rays may occur. However, extrap¬ 
olation from the LHG energy range to a higher energy 
range can be achieved by using a scaling law in the for¬ 
ward rapidity region. One possibility for such a scaling 
law is the hypothesis of limiting fragmentation |12II14| . 
which specifies that the secondary particles will approach 
a limiting distribution of rapidity in the rest frame of the 
target hadron. In this case the fragmentation of a col¬ 
liding hadron would occur independently of the center- 
of-mass energy and then the differential cross sections as 
a function of rapidity (hereafter rapidity distributions) 
in the fragmentation region, namely the forward rapidity 
region, would form a limiting distribution. 

Understanding particle production in nucleon-nucleus 
or nucleus-nucleus interactions is also of importance for 
ultrahigh-energy cosmic ray interactions, where parton 
density in nuclei is expected to be enhanced by oc 
The presence of a high gluon density in the nucleus 
is known to greatly modify the absolute yield and the 
momentum distribution of the particles that are pro¬ 
duced [T5] . 

The LHGf experiment m is designed to measure the 
hadronic production cross sections of neutral particles at 
very forward angles in p + p and proton-lead (p + Pb) 
collisions. The LHGf experiment also provides a unique 
opportunity to investigate all the effects mentioned in the 
previous paragraph, namely, the limiting fragmentation, 
the Feynman scaling im, and the high parton density 
in nuclear target. In a previous publication |18| we pre¬ 
sented the TT® production cross sections as a function of 
the transverse momentum (hereafter px distributions) in 
p + p collisions at ^/s = 7 TeV. However tests of the lim¬ 
iting fragmentation and the Feynman scaling predictions 
were not performed. Gonversely, in the analysis of this 
paper, the comparison of the LHGf data taken in p -|- p 
collisions at y/s = 2.76 and 7 TeV respectively makes it 
possible to perform these tests. In addition the analysis 
presented in this paper has updates that lead to a deeper 
understanding of forward tt® production compared to our 
previous publications |18l I19| : the upper range for px 
analysis is extended to 1.0 GeV and differential cross sec¬ 
tions as a function of longitudinal momentum (hereafter 
Pz distributions) as well as px distributions are presented. 

The paper is organized as follows. In Se c.[H| the LHGf 
detectors are described. Sections IHII and IIVI summarize 
the conditions for taking data and the MG simulation 
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methodology, respectively. In Sec. |V] the analysis frame¬ 
work and the factors that contribute to the systematic 
uncertainty of the results are explained. In Sec. |VI| the 
analysis results are presented and compared with the 
predictions of several hadronic interaction models. In 
Sec. |VII| the analysis results for p+p and p-|-Pb collisions 
are described. Finally, concluding remarks are found in 

Sec.lMIl 


II. THE LHCF DETECTOR 

Two independent detectors called LHCf Arml and 
LHCf Arm2 were assembled to study p+p and p-|-Pb col¬ 
lisions at the LHC |1D]. In p+p collisions at -^s = 7TeV, 
both LHCf Arml and LHCf Arm2 detectors were oper¬ 
ated to measure the neutral secondary particles emitted 
into the positive and negative large rapidity regions, re¬ 
spectively. In p-|-p collisions at y/s = 2.76 TeV andp-|-Pb 
collisions at y'sNN = 5.02 TeV, only the LHCf Arm2 de¬ 
tector was used to measure the neutral secondary parti¬ 
cles emitted into the negative rapidity region (the proton 
remnant side in p + Pb collisions). Here the rapidity y is 
defined as y = ta,nh~^{p^/E) [21]. 

The LHCf detectors each consist of two sampling and 
imaging calorimeters composed of 44 radiation lengths 
of tungsten and 16 sampling layers of 3 mm thick plastic 
scintillators. The transverse sizes of the calorimeters are 
20x20mm^ and 40x40mm^ for Arml, and 25x25mm^ 
and 32x32mm^ for Arm2. The smaller and larger 
calorimeters are hereafter called the Small Calorimeter 
and the Large Calorimeter, respectively. Four X-Y lay¬ 
ers of position-sensitive detectors are interleaved with 
the layers of tungsten and scintillator in order to pro¬ 
vide the transverse profiles of the showers. Scintillat¬ 
ing fiber (SciFi) belts |21] are used for Arml and sili¬ 
con microstrip sensors |23| are used for Arm2. Readout 
pitches are 1mm and 0.16 mm for Arml and Arm2, re¬ 
spectively. The Front Counters, additional components 
of the LHCf detectors, are simple thin plastic scintilla¬ 
tors (SOxSOmm^) and are installed in front of the LHCf 
calorimeters. They act as monitors for beam-beam col¬ 
lision rates with a higher detection efficiency than the 
LHCf calorimeters. 

The LHCf detectors were installed in the instrumen¬ 
tation slots of the target neutral absorbers (TANs) [21] 
located ±140 m from the ATLAS interaction point (IPl) 
in the direction of the LHCb interaction point for Arml 
and in the direction of the ALICE interaction point for 
Arm2, and at a zero-degree collision angle. The trajecto¬ 
ries of charged particles produced at IPl and directed to¬ 
wards the TANs are deflected by the inner beam separa¬ 
tion dipole magnets D1 before reaching the TANs them¬ 
selves. Consequently, only neutral particles produced at 
IPl enter the LHCf detectors. The vertical positions of 
the LHCf detectors in the TANs are manipulated so that 
the LHCf detectors cover the pseudorapidity range from 
8.4 to infinity for a beam crossing half angle of 145 prad. 


The Small Calorimeter effectively covers the zero-degree 
collision angle. Following p + Pb collision operation, the 
LHCf detectors were removed from the TAN instrumen¬ 
tation slots in April, 2013 in order to protect them from 
radiation damage when the LHC is operated at high lu¬ 
minosity. 

LHCf triggers are generated at three levels |25|. The 
first level trigger is generated from beam pickup signals 
when a bunch passes IPl. A shower trigger is then gener¬ 
ated when signals from any successive three scintillation 
layers in any calorimeter exceeded a predefined threshold. 
The shower trigger threshold is chosen to detect photons 
greater than 100 GeV with an efficiency of > 99 %. A sec¬ 
ond level trigger is generated when a shower trigger has 
occurred and the data acquisition system is activated. 
The highest level trigger, or third level trigger, is gen¬ 
erated when a specified combination of shower triggers, 
front counter triggers and data acquisition trigger has oc¬ 
curred. The live time efficiency of the data acquisition 
systems is defined as the ratio of the number of second 
level triggers to the number of shower triggers. The effi¬ 
ciency depends on the luminosity during the data taking 
period and is always less than unity due to pileup. The 
final results shown are corrected for the live time effi¬ 
ciency. 

More details on the scientific goals of the experiment 
are given in Ref. [16]. The performance of the LHCf 
detectors has been studied in previous reports pSl [25] • 


III. EXPERIMENTAL DATA TAKING 
CONDITIONS 

The experimental data used for the analysis in this pa¬ 
per were obtained at three different collision energies and 
colliding particle configurations. Data taking conditions 
are explained in the subsections below, ordered accord¬ 
ing to the dates of the operation periods with the earliest 
first. 


A. p + p collisions at •,/s = 7 TeV 

The data in p + p collisions at = 7 TeV with a zero- 
degree beam crossing angle were obtained from May 15 
to 22, 2010 (LHC Fills 1104, 1107, 1112, and 1117). The 
events that were recorded during a luminosity optimiza¬ 
tion scan and a calibration run were removed from the 
data sets for this analysis. The integrated luminosities 
for the data analysis reported in this paper were derived 
from the counting rate of the Front Counters m and 
were 2.67nb“^ (Arml) and 2.10nb“^ (Arm2) after tak¬ 
ing the live time efficiencies into account. 

Pileup interactions in the same bunch crossing may 
increase the multi-hit events that have more than one 
shower event in a single calorimeter, leading to a poten¬ 
tial bias in the momentum distributions of tt^s. The con¬ 
tamination of multi-hit events due to pileup interactions 
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is estimated to be only 0.2% and therefore produces a 
negligible effect [T^ . Detailed discussions of background 
events from collisions between the beam and residual gas 
molecules in the beam tube can be found in a previous 
report ^5] , 


B. p + Pb collisions at y^SNN = 5.02 TeV 

The data in p + Pb collisions were obtained at y^SNN = 
5.02 TeV with 145 prad beam crossing half angle and with 
only the Arm2 detector recording data on the proton 
remnant side. The beam energies were 4 TeV for protons 
and 1.58 TeV per nucleon for Pb nuclei. Because of the 
asymmetric beam energies where the proton beam travels 
at 9 = TT and the Pb beam at 9 = 0, the nucleon-nucleon 
center-of-mass in p + Pb collisions is shifted to rapidity 
—0.465 (= 1/2 X log((ApZpb)/(.^pApb)) where A and Z 
are the mass and atomic numbers, respectively |28|1. 

Data used in this analysis were taken in two differ¬ 
ent fills; during LHC Fill 3478 on January 21, 2013 and 
during LHC Fill 3481 on January 21 and 22. The in¬ 
tegrated luminosity of the data was 0.63nb“^ after cor¬ 
recting for the live time efficiencies of the data acquisition 
systems [IS]. The trigger scheme was essentially identi¬ 
cal to that used in p -l-p collisions at y/s = 7 TeV. The 
bunch spacing in p -|- Pb collisions (200 ns), which was 
smaller than the gate width for analog to digital conver¬ 
sion in the LHCf data acquisition system (500 ns) cre¬ 
ated the possibility of integrating two or at most three 
signal pulses from the pileup of successive p -1- Pb colli¬ 
sions. The probability for this to occur was estimated 
from the timing distribution for shower triggers and was 
less than 5%. Contamination by successive collisions is 
not corrected for in this study, while it is considered in 
the beam-related systematic uncertainty. The contam¬ 
ination of multi-hit events due to pileup interactions is 
negligible (0.4%). 

It should be remarked that beam divergence causes 
a smeared beam spot at the TAN, leading to a bias in 
the measured momentum distributions. The effect of a 
non-zero beam spot size at the TAN was evaluated with 
MC simulations (see Ref. jUj). This effect is taken into 
account in the final results reported for the px and pz 
distributions. 


C. p + p collisions at = 2.76 TeV 

The data in p -|-p collisions at y/s = 2.76 TeV were 
obtained with a 145prad beam crossing half angle and 
beam energy 1.38 TeV for each proton. Data used in this 
analysis were taken during LHC Fill 3563 on February 
13, 2013. The integrated luminosity for this data was 
2.36 nb“^ after correcting for the live time efficiencies of 
the data acquisition system m- The trigger scheme, 
trigger efficiency, and contamination of multi-hit events 
were mostly the same as the p -I- Pb collision data at 


y'SNN = 5.02 TeV. The effects of beam divergence were 
dealt with in the same way as was described for p -|- Pb 
collisions at ^snn = 5.02TeV (Sec. HIBl. 


IV. MONTE CARLO SIMULATIONS 
METHODOLOGY 


MC simulations have been performed in two steps: 

(I) Event generation in p -|- p and p -|- Pb collisions at 
IPl (Sec. IV AI and (H) particle transport from IPl to 
the LHCf detectors and consequent simulation of the re¬ 


sponse of the LHCf detectors (Sec. IVBl. 


MC simulation events are generated following steps (I) 
and (H) and are used for the validation of reconstruction 
algorithms, determination of cut criteria, and determina¬ 
tion of the response matrix for momentum distributions 
unfolding. Conversely, MC simulations that are used only 
for comparison with the measurement results in Sec. |VI| 
are limited to step (I) only, since the final px and p^ dis¬ 
tributions in Sec. |VI| are already corrected for detector 
response and eventual reconstruction bias. The statisti¬ 
cal uncertainties of the MC simulations used in this paper 
are negligibly small compared to the statistical uncertain¬ 
ties of the LHCf data. 


A. Collision event modeling 

Collision event modeling oi p + p hadronic interactions 
at y/s = 2.76 and 7TeV are simulated and the resulting 
fluxes of secondary particles are generated with several 
event generators: dpmjet 3.06 m, QGSJET H-04 122], 
SIBYLL 2.1 m , EPOS LHC m and PYTHIA 8.185 [251123]. 
Hereafter the version number for these event generators 
is omitted for simplicity, unless otherwise noted. 

In the analysis of this paper, we use the integrated 
interface CRMC 1.5.3 m for executing the first four event 
generators, whereas the fifth event generator, pythia, 
serves as its own front end for the generation of proton- 
proton hadronic interaction events. 

Events in p -|- Pb collisions are divided into two cat¬ 
egories according to the value of the impact parameter: 
(1) general hadronic interactions and (2) ultraperipheral 
collisions (UPCs). Category (1) occurs when the im¬ 
pact parameter between p and Pb is smaller than the 
sum of their radii. These inelastic p -|- Pb interactions 
at .^snn = 5.02 TeV are simulated using the hadronic 
interaction models dpmjet, QGSJET, and EPOS with the 
CRMC interface. SIBYLL was not used because it only 
supports nuclei lighter than Fe. pythia also does not 
support heavy ion collisions and thus was also not used 
for p -I- Pb collisions. 

Category (2) p -|- Pb UPCs occur when the impact 
parameter is larger than the sum of p and Pb radii. 
The UPC events are simulated by the combination of 
STARLIGHT [25] for the virtual photon flux, SOPHIA 
2.1 [39] for low-energy photon-proton interactions, and 
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either dpmjet 3.05 m or PYTHIA 6.428 |3S] for high- 
energy photon-proton interactions. The UPC simula¬ 
tion distributions used in this analysis are taken from 
the average of two UPC simulations; one using dpmjet 
3.05 and the second using PYTHIA 6.428 for the high- 
energy photon-proton interaction. Differences between 
these two UPC simulations are taken into account as 
a systematic uncertainty in the UPC simulation. See 
Ref. 130] for more details. 

In both p+p and p-|-Pb collisions, the MC events used 
for the determination of the response matrix for unfolding 
the momentum distributions (Sec. VBl are simulated by 
PYTHIA at the requisite beam energies. A single 7r° with 
energy larger than 100 GeV and possible associated back¬ 
ground particles are selected from the secondary particles 
produced. There is no significant dependence of the un¬ 
folding performance on the choice of event generator for 
the MC simulation events that are used for the response 
matrix. This was verified by repeating event simulations 
with three of the event generators; dpmjet, PYTHIA and 
EPOS. 


In all of the MC simulations, the 7r°s from short-lived 
particles that decay within Im of IPl, e.g. rj, p, ui, etc. 
(< 10% for each relative to all 7r°s), are accounted for 
consistently in the treatment of LHCf data. The 145 prad 
beam crossing half angle is also taken into account for 
p+p collisions at y/s = 2.76 TeV and for p-|-Pb collisions 
at Usnn = 5.02 TeV. 


B. Simulation of particle transport from IPl to the 
LHCf detector and of the detector response 


Transport of secondary particles inside the beam pipe 
from IPl to the TAN, the electromagnetic and hadronic 
showers produced in the LHCf detector by the particles 
arriving at the TAN and the detector response are simu¬ 
lated with the COSMOS and EPICS libraries El- 

Secondary particles produced by the interaction be¬ 
tween IPl collision products and the beam pipe are also 
taken into account in this step. The secondary particles 
from beam pipe interaction events generally have energy 
well below 100 GeV and thus provide no bias to the mo¬ 
mentum distributions of collision events that focus only 
on energies above 100 GeV. The survey data for detector 
position and random fluctuations due to electrical noise 
are also taken into account in this step. See Ref. [TB] for 
more details. 


V. ANALYSIS FRAMEWORK 

A. event reconstruction and selection 

The standard reconstruction algorithms consist of four 
steps: hit position reconstruction, energy reconstruction, 
particle identification, and 7r° event selection. 


1. Position reconstruction 


Hit position reconstruction starts with a search for 
multi-hit and single hit events. A multi-hit event is de¬ 
fined to have more than one photon registered in a single 
calorimeter. A single-hit event is defined to have a single 
hit in each of the two calorimeters in a given detector, 
Arml or Arm2. 

Therefore multi-hit event candidates should have two 
or more distinct peaks in the lateral-shower-impact- 
distribution of a given calorimeter and are then iden¬ 
tified using the TSpectrum algorithm |42| implemented 
in ROOT |33|. TSpectrum provided the basic function¬ 
ality for peak-finding in a spectrum with a continuous 
background and statistical fluctuations. 

The MC simulation estimated efficiencies for identi¬ 
fying multi-hit events are larger than 70% and 90% 
for Arml and Arm2, respectively [3S|. Given the list 
of shower peak position candidates that have been ob¬ 
tained above, the lateral distributions are fit to a Loren- 
zian function j44] to obtain more precise estimates of the 
shower peak positions, heights, and widths. In the case 
of multi-hit events, two peaks are fit using superimposed 
Lorenzian functions. Multi-hit events with three or more 
peaks are rejected from the analysis. Conversely, single¬ 
hit events, not having two or more identifiable peaks 
in a single calorimeter but having a single hit in each 
calorimeter are correctly selected with an efficiency bet¬ 
ter than 98% for true single-photon events with energy 
greater than 100 GeV for both Arml and Arm2. 


2. Energy reeonstruetion 


The photon energy is reconstructed using the measured 
energy deposited in the LHCf calorimeters. The charge 
information in each scintillation layer is first converted to 
a deposited energy by using the calibration coefficients 
obtained from the SPS electron test beam data taken 
below 200 GeV |26]. The sum of the energy deposited 
in the 2nd to 13th scintillation layers is then converted 
to the primary photon energy using an empirical func¬ 
tion. The coefficients of the function are determined from 
the response of the calorimeters to single photons using 
MC simulations. Corrections for shower leakage effects 
and the light-yield collection efficiency of the scintillation 
layers are carried out during the energy reconstruction 
process |20|. In the case of multi-hit events, the recon¬ 
structed energy based on the measured energy deposited 
is split into two energies, primary and secondary. Frac¬ 
tions of the energy for the primary and secondary hits 
are determined according to the peak height and width 
of the corresponding distinct peaks in the lateral-shower- 
impact-distribution . 
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3. Particle identification 

Particle identification (PID) is applied in order to effi¬ 
ciently select pure electromagnetic showers and to reduce 
hadron (predominantly neutron) contamination. PID in 
the study of this paper depends only on the parameter 
Lgo%. ^ 90 % is defined as the longitudinal distance, in 
units of radiation length (Xg), measured from the 1st 
tungsten layer of the calorimeter to the position where 
the energy deposition integral reaches 90% of the total 
shower energy deposition. Events with an electromag¬ 
netic shower generally have a Lgo% value smaller than 
20 Xq, while events with a hadronic shower generally have 
LgQ% larger than 20 Xq. The threshold LgQ% value as a 
function of the photon energy is defined in order to keep 
the 7r° selection efficiency at 90 % over the entire energy 
range of the individual photons. PID criteria are deter¬ 
mined by MC simulations for each calorimeter. 

4- 7r° event selection 

The 7r° are then identified by their decay into two pho¬ 
tons, leading to the distinct peak in the invariant mass 
distribution around the rest mass. The invariant mass 
of the two photons is calculated using the reconstructed 
photon energies and incident positions. The 7r° events 
used in the analysis of this paper are classified into two 
categories: Type-I 7r° and Type-II 7r° events. A Type- 
I event is defined as having a single photon in each of 
the two calorimeters of Arml or Arm2 (the left panel of 
Fig. 0 . A Type-II event is defined as having two photons 
in the same calorimeter (the right panel of Fig. 0. Note 
that Type-II events were not used in the previous analy¬ 
ses HHmii, and thus are taken into account for the first 
time in this paper. As detailed in Sec. |V B| the phase 
spaces covered by Type-I and Type-II events are comple¬ 
mentary. In particular, the inclusion of Type-II events 
extends the px upper limit for analysis from 0.6 GeV in 
the previous analyses to 1.0 GeV. 

Figure 0 shows the reconstructed two-photon invari¬ 
ant mass distributions of LHGf data in the rapid¬ 

ity range 8.8 < y < 10.8. The left and right panels of 
Fig.0 show the distributions for Type-II events in the 
Arm2 small calorimeter and Arm2 large calorimeter re¬ 
spectively. The sharp peaks around 135 MeV are due to 
7r° events. The distributions in Fig. 0are based only on 
data from p + p collisions at ^/s = 7 TeV during LHG Fill 
1104. Similar invariant mass distributions are obtained 
from other fills and from Arml. Kinematic quantities of 
the TT^s (four-momenta, pxj Pz and rapidity) are recon¬ 
structed by using the photon energies and incident posi¬ 
tions measured by the LHGf calorimeters, and are used 
for producing the px and p^ distributions. The projected 
position of the proton beam axis on the LHGf detector 
(beam center) is used in order to derive the correct px 
and Pz values of each event. The beam center position 
is obtained from the LHGf position-sensitive detectors of 



FIG. 1. (color online). Observation of 7r° decay by a LHGf 
detector. Left: Type-I 7r° event having one photon entering 
each calorimeter. Right: Type-II tt® event having two photons 
entering one calorimeter, here entering the small calorimeter. 


Arml and Arm2 for each fill. 




FIG. 2. (color online). Reconstructed invariant mass dis¬ 
tributions in p-I-p collisions at -^5 = 7 TeV. Left: Type-II 
7r° events in the Arm2 small calorimeter. Right: Type-II tt® 
events in the Arm2 large calorimeter. The solid curves show 
the best-fit composite physics model to the invariant mass 
distributions. 

The TT^ event selection criteria that are applied prior to 
the reconstruction of the ir^ kinematics are summarized 
in Table. Type-I events accompanied by at least one 
additional background particle in one of the two calorime¬ 
ters (usually a photon or a neutron) and not originating 
in a 7r° decay are denoted as multi-hit 7r° events and are 
rejected as background events. Similarly, Type-II events 
accompanied by at least one additional background par¬ 
ticle in the calorimeter used for tt^ identification are re¬ 
jected. Figure 0 shows diagrams of all types of multi-hit 
events that are rejected. Panels (a) and (b) show the 
multi-hit Type-I 7r° events and panels (c) and (d) show 
the multi-hit Type-II 7r° events. Red and green arrows 
indicate a background particle not originating in a 
decay and two photons originating in a decay, respec¬ 
tively. The final inclusive production rates reported in 
this paper are corrected for these cut efficiencies and will 
be discussed in Sec. EB 
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FIG. 3. (color online). Diagrams of all multi-hit events that 
are rejected. Panels (a) and (b) show the multi-hit Type-I 7r° 
events and panels (c) and (d) show the multi-hit Type-II tt® 
events. Red and green arrows indicate a background particle 
not originating in a tt® decay and two photons originating in 
a 7r° decay, respectively. 


is performed over the two photon invariant mass range 
0.08 < < 0.18 GeV. The 7r° signal window is defined 

by the two dashed vertical lines in Fig. that are placed 
±3 ct from the mean value. Here the mean value and the 
standard deviation are obtained from the best-fit asym¬ 
metric Gaussian distribution. The background window is 
defined as the region within ±6cr distance from the peak 
value and excluding the tt® signal window. The fraction 
of the background component included in the 7r° signal 
window can be estimated using the ratio of the integral 
of the best-fit third order Ghebyshev function over the 7r° 
signal window divided by the integral over the 7r° signal 
and background windows. The width of the asymmet¬ 
ric Gaussian function comes from the detector response, 
predominantly from shower leakage near the edges of the 
calorimeters. The reconstructed energy is corrected for 
shower leakage. 


TABLE I. Summary of criteria for selection of the 7r° sample. 
Type-I 7r° events 

Incident position Within 2 mm from the edge of calorimeter 
Energy threshold ifphoton >100 GeV 
Number of hits Single-hit in each calorimeter 
PID Photon-like in each calorimeter 

Type-II 7r° events 

Incident position Within 2 mm from the edge of calorimeter 
Energy threshold ifphoton >100 GeV 
Number of hits Two hits 
PID Photon-like 


B. Corrections for experimental effects 

The raw px distributions of tt^s are corrected 

for: (1) contamination by background events, (2) recon¬ 
struction inefficiency and the smearing caused by finite 
position and energy resolutions, (3) geometrical accep¬ 
tance and the branching ratio of 7r° decay, and (4) the 
efficiency of the multi-hit 7r° cut. We now discuss each 
of these corrections in some detail. 


1. Background contamination 

First, the background contamination of the events 
from hadronic events, and from the coincidence of two 
photons not originating from the decay of a single 7r° are 
estimated using a sideband method [18]. As shown in 
Fig.jg for instance, the reconstructed two-photon invari¬ 
ant mass distributions of LHGf data are fit to a com¬ 
posite physics model (solid blue curve). The model con¬ 
sists of an asymmetric Gaussian distribution for the 7r° 
signal component and a third order Ghebyshev polyno¬ 
mial function for the background component. The fit 


2. Reconstruction inefficiency and smearing in position and 
energy resolution 

Second, a spectrum unfolding is performed to simul¬ 
taneously correct for both the reconstruction inefficiency 
and the smearing caused by finite position and energy res¬ 
olution. The contamination by background events that 
has been estimated by the sideband method is taken 
into account in the unfolding process. We follow basi¬ 
cally same unfolding procedure as in the previous analy¬ 
ses nmn], although the unfolding algorithm is based on 
a fully Bayesian unfolding method |45| instead of an iter¬ 
ative Bayesian unfolding method |3B|. The calculation of 
the “a posteriori” probability in multi-dimensional space 
(the measured spectrum multiplied by the true spec¬ 
trum) is achieved using a Markov Ghain Monte Garlo 
simulation m- The convergence of the Markov Ghain 
Monte Garlo simulation is ensured by the Gelman-Rubin 
test jlH]. Production of the MG events used for the cal¬ 
culation of the response matrix for the unfolding is ex¬ 
plained in Sec. |IV A[ 

3. Geometric acceptance and branching ratio corrections 

Thirdly, the limiting aperture of the LHGf calorimeters 
is estimated by using MG simulations. The procedure for 
performing MG simulations is given in Ref. m- Figured 
shows the acceptance efficiency as a function of the 
Pz and px- The acceptance efficiency has been obtained 
by taking the ratio of the p^rPn distribution of tt^s that 
are within the aperture of the LHGf calorimeters divided 
by the distribution of all simulated tt^s. The fiducial 
cuts [TS] and reconstructed energy cut (both of the 7r° 
decay photons must have E > 100 GeV) are also applied 
to the accepted 7r° events. Dashed curves in Fig. [^in¬ 
dicate lines of constant 7r° rapidity. The acceptance effi¬ 
ciencies in Fig. [^ are purely kinematic and do not depend 
upon a particular hadronic interaction model. The aper- 






















ture correction is achieved by dividing, point by point, 
the distributions before the acceptance correction by the 
acceptance efficiency. The branching ratio inefficiency is 
due to 7r° decay into channels other than two photons. 
The branching ratio for 7r° decay into two photons is 
98.8% and is taken into account by increasing the 7r° 
acceptance efficiency by 1.2%. 


4-. Loss of events due to the multi-hit tt® cut 


Fourth, in order that the reported 7r° distributions 
represent inclusive cross sections it is necessary to cor¬ 
rect the data for the loss of events due to the multi¬ 
hit cut (Sec. |V A4|). The correction factor is defined 

as /multihit ^ ^single ^ ^^single ^ 

and are the number of expected multihit and 

single-hit 7r° events in the Tth bin respectively. The 
factors estimated using hadronic interac¬ 

tion models introduced in Sec. |IV A] and are in the range 
1.0 < ^2.1 over all the px and bins. LHCf px 

and pz distributions are then multiplied by the average 
of these factors for the various interaction models and 
their contribution to the systematic uncertainty is de¬ 
rived from the observed variations amongst the interac¬ 
tion models. Consequently, the single-hit 7r° distributions 
are corrected to represent inclusive 7r° production distri¬ 
butions. All the procedures just described have been ver¬ 
ified using the MC simulations introduced in Sec. IV A 



1000 2000 3000 4000 1000 2000 3000 4000 


P, [GeV] [GeV] 

FIG. 4. (color online). The acceptance map of detection by 
the LHCf detectors in Pz-px phase space: Arml Type-I (left 
top), Arml Type-II (right top), Arm2 Type-I (left bottom), 
and Arm2 Type-II (right bottom). The fiducial area cuts and 
energy threshold (Aphoton > 100 GeV) are taken into account. 
Dashed curves indicate lines of constant rapidity tt^s, y = 8.8, 
9.0, and 10.0 reading from top to bottom. 


C. Systematic uncertainties 

Systematic uncertainties are determined by three fac¬ 
tors: (1) possible biases in event reconstruction, (2) un¬ 
certainty of the LHC machine conditions, and (3) an in¬ 
teraction model dependence. 

1. Systematic uncertainties in event reconstructions and 
unfolding of distributions 

Uncertainties related to biases in event reconstruction 
are mainly due to five causes: (1) single-hit/multi-hit 
separation, (2) PID, (3) energy scale uncertainty, (4) 
position-dependent corrections for both shower leakage 
and the light yield of the calorimeters, and (5) unfolding 
of distributions. For the first four terms, we follow the 
same approaches to estimate the systematic uncertainties 
as we used in the previous study |18| . 

Concerning the unfolding process, the uncertainty is 
estimated by adding the following three components in 
quadrature. First, the uncertainty due to a possible de¬ 
pendence of the unfolding procedure on the shape of 
the px or distributions to be unfolded is estimated 
from MC simulations; we estimate the variation of the 
ratios of the unfolded distributions to the true distribu¬ 
tions among the three true distributions predictions by 


DPMJET, QGSJET, and EPOS. The second component is 
a dependence of the unfolding procedure on the event 
generator used in the generation of the response ma¬ 
trix for unfolding, which is negligible as we mentioned 
in Sec. |IV A| Finally, the third component is the system¬ 
atic uncertainty in the unfolding algorithm itself. This 
is evaluated by comparing two unfolded distributions, 
one obtained by a fully Bayesian unfolding method and 
the second obtained by the iterative Bayesian unfolding 
method. The uncertainty in the first component is 10 % 
over the all px and pz bins, and the uncertainties in the 
other two components make no significant contribution. 
Thus we assign 10 % for the systematic uncertainty in the 
unfolding of px and pz distributions. 


2. Systematic uncertainties in the LHC machine conditions 

The LHC machine conditions introduce systematic un¬ 
certainties in beam position and luminosity. The beam 
position at the LHCf detectors varies from fill to fill ow¬ 
ing to variations of the beam transverse position and 
the crossing angles at IPl. The beam center positions 
at the LHCf detectors obtained for LHC Fills 1089 to 
1134 by the LHCf position-sensitive detectors and by the 
beam position monitors (BPMSW) installed ±21 m from 
IPl jini are consistent with each other within ±lmm. 
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The systematic shifts to the px and pz distributions are 
then evaluated by taking the ratios of distributions with 
the beam center displaced by ±lmm to distributions 
with no displacement present. The evaluated systematic 
shifts to the pt and pz distributions are 5-20 % depend¬ 
ing on the Pt and pz values. 

The uncertainty in the luminosity depends on the col¬ 
lision configuration. For the data in p + p collisions at 
y/s = 7TeV, the luminosity value used for the analysis 
is derived from the counting rate of the Front Counters. 
Considering the uncertainties in both the calibration of 
the Front Counters ±3.4% and in the beam intensity 
measurement ±5.0% during the Van der Meer scans, we 
estimate an uncertainty of ±6.1% in the luminosity for 
p + p collisions at ^/s = 7 TeV m For the p + p col¬ 
lision data at ^/s = 2.76 TeV and p ± Pb collision data 
at y'sNN = 5.02 TeV, LHCf data were taken simultane¬ 
ously with data taken by the ATLAS experiment. The 
luminosity values used for this data analysis were then 
provided by the LHCf Front Counters and also by the 
ATLAS collaboration. 

The luminosity uncertainties in p+p collisions at ^/s = 
2.76 TeV and in p±Pb collisions at = 5.02 TeV are 
estimated to be ±3.1 % [M] and ±20% |in|, respectively. 

Pileup of successive p ± Pb collisions due to the small 
bunch spacing (200 ns) relative to the data acquisition 
time (500ns) amounts to < 5% systematic uncertainty of 
Pt and pz distributions (see Sec. IIIB I, and may provide 
a slight shift of the absolute normalization for the px and 
Pz distributions. This effect is not corrected for in this 
study, but is taken into account as uncertainty related to 
the LHC machine condition. 


TABLE II. Summary of the systematic uncertainties. Numer¬ 
ical values indicate the maximum variation of bin contents 
in the px and Pz distributions due to systematic uncertain¬ 
ties. Note that the uncertainty in contamination of successive 
p -I- Pb collisions and in UPC 7r° simulation pertain only to 
p + Pb collisions. 


Single/multihit selection 

±3% 

Particle identification 

±(0-20%) 

Energy scale 

±(5-20%) 

Position-dependent correction 

±(5-30%) 

Unfolding 

±(5-10%) 

Offset of beam axis 

±(5-20%) 

Luminosity (p -|- p at 7 TeV) 

±6.1% 

Luminosity (p-|-p at 2.76 TeV) 

±3.1% 

Luminosity (p -|- Pb at 5.02 TeV) 

±20% 

Contamination of successive p -|- Pb collisions 

<5% 

Multihit 7r° correction 

<10% 

UPC 7r° simulation 

±(3-20%) 


the 1st paragraph of Sec. |V C 2[ 5-20 % shifts in px or pz 
distributions are obtained, (7) Luminosity uncertainty is 
explained in the 2nd paragraph of Sec. V C 2 (8) Con¬ 
tamination of successive p ± Pb collisions is explained 
in the 3rd paragraph of Sec. |V C 2[ (This uncertainty 
is due to contamination and thus only a positive error 
is quoted.) (9) The uncertainty in multihit 7r° events 
±10%, and (10) the uncertainty in UPC ±(3-20%) are 
found in Sec. |V C 3| Table summarizes the systematic 
uncertainties of the tt^ px and pz distributions. 


VI. ANALYSIS RESULTS 


3. Systematic uncertainties depending on the interaction 
models used in the MC simulations 

The analysis in this paper unavoidably relies on the 
predictions given by MC simulations. First, we correc t 
LHCf data for the loss of multi-hit 7r° events (Sec. V B 41. 
The correction factors g. systematic un¬ 

certainty of less than 10 % among the hadronic interac¬ 
tion models. Second, for p ± Pb collisions only, the con¬ 
tamination from UPC induced 7r° events in LHCf data 
is derived from MC simulations (Sec. |I V A I. The com¬ 
parison of the predicted px and pz distributions of 7r°s 
between two UPC MC simulations, one using dpmjet 
3.05 and the other one using pythia 6.428 for high- 
energy photon-proton interaction, show a systematic un¬ 
certainty of roughly 3-20%. 

In summary, there are 10 systematic uncertainties. 
The first four (1) Single/multihit selection, (2) PID, (3) 
energy scale and (4) position-dependent correction are 
explained in Ref. m and we follow the same approaches 
as we used in Ref. m- The remaining six systematic 
uncertainties and the text containing their explanations 
are: (5) Unfolding uncertainty is explained and evalu¬ 
ated in Sec. V C 1 (6) Offset of beam axis is explained in 


A. Results in p-|-p collisions at y's = 7 TeV 

The inclusive production rate of neutral pions as a 
function of px and Pz is given by the expression m 


1 „d^(J 1 (fN{p-r,y) 


1 ^d‘^N{pT,Pz) 


Ninel 2TrpTdpTdy Ninel 2npTdpTdp, 


CTinel ^ dp^ 

( 1 ) 

CTinei is the inelastic cross section for p ± p collisions at 
= 7 TeV. Ed^a/dp^ is the inclusive cross section for 
TT^ production. The number of inelastic collisions, Mnei, 
used for the production rate normalization is calculated 
from Ninel = CTinei / dldt, taking the inelastic cross section 
o'inei = 73.6 mb m- The uncertainty in o’inei is estimated 
to be ±3.0 mb by comparing the values of o’inei reported 
in Refs. [5DH5?| . 

Us ing th e integrated luminosities f Cdt, reported in 
Sec. HI A Ainei is (2.67 ± 0.11) x 10® for Arml and 
(2.10 ±0.09) X 10® for Arm2. d'^N{p^,y) is the num¬ 
ber of TT^s produced within the transverse momentum 
interval dpx and the rapidity interval dy. Similarly 
d^N{pT,Pz) is the number of 7r®s produced within dpx 
and the longitudinal momentum interval dpz- 

Experimental px and pz distributions measured inde- 
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pendently with the Arml and Arm2 detectors are com¬ 
bined following a pull method |54| and the final px and pz 
distributions are then obtained by minimizing the value 
of the chi-square function defined by 


X 


= EE 


RTia + Sa,^) - R\ 


comb 


Xpenalty (^) 


where the index i represents the px or pz bin number run¬ 
ning from 1 to n (the total number of px or pz bins), and 
the index a indicates the type of distributions: a = 1 
Arml Type-I events, a = 2 Arml Type-II events with 
the Large Calorimeter, a = 3 Arm2 Type-I events, a = 4 
Arm2 Type-II events with the Small Calorimeter, and 
0 = 5 Arm2 Type-II events with the Large Calorimeter. 
Note that Arml Type-II events with the Small Calorime¬ 
ter are not used for this analysis since the energy recon¬ 
struction accuracy for these events is still being investi¬ 
gated. is the inclusive production rate in the Ah bin 
of the oth distributiom which corresponds to the second 
and third terms Eq. Q. is the inclusive produc¬ 

tion rate in the Ah bin obtained by combining all R‘^l’s 
for a =1-5. CTa.i is the uncertainty of The aa,i 

are calculated by quadratically adding the statistical un¬ 
certainty and the systematic uncertainty in the energy 
scale. The energy scale uncertainty has been estimated 
with test beam data taken at the SPS and is uncorrelated 
bin-by-bin unlike the other systematic uncertainties |18| . 
The systematic correction Sa^i modifies the number of 
events in the Ah bin of the ath distribution: 




i=i 


■ 

J a,I 


(3) 


The coefficient ^ is the systematic shift of the Ah bin 
content of the ath distribution due to the jth systematic 
uncertainty term. The systematic uncertainty consists of 
seven uncertainties related to the single-hit/multihit sep¬ 
aration, the PID, the energy scale (owing to the invari¬ 
ant mass shift of the measured 7r° events), the position- 
dependent correction, the unfolding procedure, the beam 
center position, and the loss of multihit 7r° events. These 
uncertainties are assumed to be fully uncorrelated be¬ 
tween the Arml and Arm2 detectors, while correlations 
between Type-I and Type-II events and bin-bin corre¬ 
lations have been accounted for. The coefficients e^, 
which should follow a Gaussian distribution, can be var¬ 
ied, within the constraints of the penalty term given by 

Xpenalty EE kilt (4) 

j=l a=l 

to achieve the minimum value for each chi-square test. 
Note that the uncertainty in the luminosity determina¬ 
tion, ±3.1 %-±20%, not included in Eq. ^ and Eq. Q, 
can cause independent shifts of all the px and pz distri¬ 
butions. 


The LHCf px distributions (filled circles) are obtained 
from the best-fit and are shown in Fig.[^ The 68 % 

confidence intervals incorporating the statistical and sys¬ 
tematic uncertainties, except for the luminosity uncer¬ 
tainty, are indicated by the error bars. LHCf px dis¬ 
tributions are corrected for the influences of the detec¬ 
tor response, event selection efficiencies and geometrical 
acceptance efficiencies, and thus LHCf px distributions 
can be compared directly to the predicted px distribu¬ 
tions from hadronic interaction models. For comparison, 
the predictions from various hadronic interaction mod¬ 
els are also shown in Fig. dpmjet (solid red line), 
QGSJET (dashed blue line), SIBYLL (dotted green line), 
EPOS (dashed-dotted magenta line), and PYTHIA (default 
parameter set, dashed-double-dotted brown line). For 
these hadronic interaction models, the inelastic cross sec¬ 
tion used for the production rate normalization is taken 
from the predefined value in each model. 

Figure presents the ratios of the inclusive produc¬ 
tion rates predicted by the hadronic interaction models 
listed above to those obtained by LHCf data. Shaded ar¬ 
eas have been taken from the statistical and systematic 
uncertainties. In Fig. the denominator and the nu¬ 
merators, namely the inclusive production rate for LHCf 
data and for the hadronic interaction models, respec¬ 
tively, are properly normalized by the inelastic cross sec¬ 
tion for each, thus we do not apply any other normaliza¬ 
tion to the ratios. The inclusive production rates of 7r°s 
measured by LHCf and the ratios of 7r° production rate 
of MC simulation to data are summarized in Appendix. 

In the comparisons in Fig. and QGSJET has good 
overall agreement with LHCf data, while EPOS produces 
a slightly harder distribution than the LHCf data for 
Px > 0.5 GeV. These two models are based on the 
parton-based Gribov-Regge approach and are 

tuned by using the present LHC data (ALICE, ATLAS, 
CMS, and TOTEM) IM] . The prediction of SIBYLL 
agrees well with the LHCf data for 8.8 < y < 9.2 and 
Px < 0.4 GeV, while the absolute yield of SIBYLL is about 
half that of the LHCf data for y > 9.2. The predic¬ 
tions of DPMJET and PYTHIA are compatible with LHCf 
data for 9.0 < y < 9.8 and px < 0.2 GeV, while for 
Px > 0.2 GeV they become significantly harder than both 
LHCf data and the other model predictions. Generally 
the harder distributions appearing in SIBYLL, dpmjet, 
and PYTHIA can be attributed to the baryon/meson pro¬ 
duction mechanism that is used by these models. For 
example the popcorn approach |5ll|58] implemented in 
the Lund model is known to produce hard distributions 
of forward mesons |59| . Indeed, by only changing the tun¬ 
ing parameters of the popcorn approach in dpmjet one 
obtains softer meson distributions and consequently px 
distributions that are compatible with LHCf data. How¬ 
ever such a crude tune may bring disagreements between 
the model predictions and other experimental results, e.g. 
forward neutron pz and px distributions. 

The LHCf Pz distributions are shown in Fig. The 
Pz distributions predicted by various hadronic interaction 
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FIG. 5. (color online). LHCf px distributions (filled circles) in p-|-p collisions at y's = 7TeV. Error bars indicate the total 
statistical and systematic uncertainties. The predictions of hadronic interaction models are shown for comparison: dpmjet 
(solid red line), QGSJET (dashed blue line), sibyll (dotted green line), epos (dashed-dotted magenta line), and pythia (dashed- 
double-dotted brown line). 


models are also shown in Fig. Figure presents the 
ratios of Pz distributions predicted by the hadronic in¬ 
teraction models to the LHCf Pz distributions. Shaded 
areas have been taken from the statistical and systematic 
uncertainties. The same conclusions for the comparisons 
are obtained as those found for Fig. and There is 
again an overall agreement between LHCf data and the 
QGSJET prediction, especially for 0.0 < pt < 0.2 GeV. 
The EPOS prediction is compatible with LHCf data for 
Pt < 2 TeV, while showing a hard slope for px > 2 TeV 
in allpT regions. The predictions by dpmjet and pythia 
agree with LHCf data for px < 0.2 GeV and pz < 1.6 TeV, 
while showing a harder distribution for the higher pz re¬ 
gions. SIBYLL predicts a smaller production of 7r°s for 
px < 0.2 GeV and becomes similar with dpmjet and 
PYTHIA with increasing px. 


B. Results in p + p collisions at = 2.76 TeV 

The inclusive production rates of 7r°s as a function of 
Px and Pz are given by Eq. Q. Using the inelastic cross 
section tTinei = (62.5 ± 5.0) mb and the integrated 
luminosities reported in Sec. [Ill C| Ni^ei is calculated as 
(1.60 ± 0.13) X 10®. The uncertainty on dinei is estimated 
by comparing the tTinei value with the present experimen¬ 
tal result m- Note that only the LHCf Arm2 detector 
was operated in p+p collisions at y/s = 2.76 TeV and that 
only Type-I events are used for the analysis since Type-II 
event kinematics are outside the calorimeter acceptance 
for 75 = 2.76 TeV. 

LHCf Px distributions are shown in Fig. The px dis¬ 
tributions predictions for the hadronic interaction mod¬ 
els are also shown in Fig. for comparison. Figure 
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FIG. 6. (color online). Ratios of LHCf px distributions to the px distributions predicted by hadronic interaction models in 
p + p collisions at ^/s = 7TeV are shown by solid red line (dpmjex), dashed blue line (qgsjex), dotted green line (sibyll), 
dashed-dotted magenta line (epos), and dashed-double-dotted brown line (pyxhia). Shaded areas indicate the range of total 
uncertainties of the LHCf px distributions. 


presents the ratios of px distributions predicted by the 
hadronic interaction models to the LHCf px distribu¬ 
tions. QGSJET provides the best agreement with LHCf 
data, although it is slightly softer than the LHCf data 
for y > 9.2. The prediction of EPOS shows a harder be¬ 
havior than both QGSJET and LHCf data, sibyll tends 
to have generally a smaller yield and a harder dis¬ 
tribution compared to QGSJET and EPOS, leading to the 
smaller and larger yields with respect to LHCf data in 
the Px regions below and above 0.1 GeV. dpmjet and 
PYTHIA predict larger 7r° yields than both LHCf data 
and other models over the entire rapidity range. The 
same discussion on the popcorn model in the previous 
Section [Vl A| can be applied to the predictions of sibyll, 
DPMJET, and PYTHIA. 


LHCf Pz distributions are shown in Fig. Figure [l^ 

presents the ratios of p^ distributions predicted by the 
hadronic interaction models to LHCf Pz distributions. 
The same tendencies found in Fig.|^are present here also 
QGSJET gives the best agreement for 0.0 < px < 0.4 GeV 
and EPOS has a harder behavior especially for 0.2 < px < 
0.4 GeV. The predictions of dpmjet and PYTHIA are sig¬ 
nificantly harder than LHCf data for px < 0.4 GeV and 
show poor overall agreement with LHCf data. This can 
be explained by the popcorn model in a way similar to 
the harder px distributions of the SIBYLL, dpmjet and 
PYTHIA models found in Fig.j^and the preceding Section. 
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FIG. 7. (color online). LHCf Pz distributions (filled circles) in p-|-p collisions at yT = 7TeV. Error bars indicate the total 
statistical and systematic uncertainties. The predictions of hadronic interaction models are shown for comparison: dpmjet 
(solid red line), QGSJET (dashed blue line), sibyll (dotted green line), epos (dashed-dotted magenta line), and pythia (dashed- 
double-dotted brown line). 


C. Results in p -f Pb collisions at y^SNN = 5.02 TeV 

The inclusive tt® production rate in p + Pb collisions is 
given as 


1 ^ 1 d^iVPP'^(pT,yiab) 

^P^Pb ^p3 j^pPh 2TTpTdpTdyuh 

^ 1 ^ d^iVPP^pT,Pz) . 5 , 

‘^irpTdpTdp^ 

where is the inelastic cross section, Ed^a^^^/dp^ is 
the inclusive cross section of production in p -|- Pb col¬ 
lisions at y^NN = 5.02 TeV, and r/iab is the rapidity in the 
detector reference frame. The number of inelastic p-|-Pb 
collisions, used for normalizing the production 

rates is calculated from / Cdt, assuming the 

inelastic p-|-Pb cross section = (2.11 ± 0.11) b [^ . 
The value for crPjgJ’ is derived from the inelastic p + p 
cross section and the Glauber multiple collision 

model EU Eg. The uncertainty on is estimated 

by comparing the value with other calculations and 
experimental results presented in Refs. |631 [51] . Using 


III 


<1 


the integrated luminosities described in Sec. 
is (9.33 ± 0.47) x 10^. Note that only the LHCf Arm2 
detector (proton remnant side) was operated in p -|- Pb 
collisions at ^snn = 5.02 TeV. 

Figure shows LHCf px distributions with both 
statistical and systematic errors (filled circles and er¬ 
ror bars). The px distributions in p -|- Pb collisions at 
i^snn = 5.02 TeV predicted by the hadronic interaction 
models, dpmjet (solid red line), QGSJET (dashed blue 
line), and EPOS (dotted magenta line), are also shown in 
the same figure for comparison. The expected UPC con¬ 
tribution discussed in Sec. II V Al is added to the hadronic 
interaction model predictions for consistency with the 
treatment of LHCf data, and the UPC px distribution 
is shown for reference (dashed-double-dotted green line). 

In Fig. DPMJET shows good agreement with LHCf 
data at — 8.8 > piab > —10.0 and px < 0.3 GeV, while 
showing a harder behavior at — 8.8 > j/iab > —9.2 and 
Px > 0.5 GeV. QGSJET and epos predict relatively sim¬ 
ilar distributions to each other and show better agree¬ 
ment with LHCf data for px > 0.4 GeV than dpmjet. 
The characteristic bump at piab > — 9.8 and 0.1 < px ^ 
0.2 GeV, which is absent in p-|-p collisions, originates 
from the channel 7 -I- p —>■ + p via baryon resonances 
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FIG. 8. (color online). Ratios of LHCf Pz distributions to the pz distributions predicted by hadronic interaction models in 
p + p collisions at ^ = 7TeV are shown by solid red line (dpmjet), dashed blue line (qgsjet), dotted green line (sibyll), 
dashed-dotted magenta line (epos), and dashed-double-dotted brown line (pythia). Shaded areas indicate the range of total 
uncertainties of the LHCf Pz distributions. 


in UPCs. In fact the UPC simulation reproduces such 
a bump. Figure presents the ratios of LHCf px dis¬ 
tributions to the pt distributions predicted by hadronic 
interaction models taking the UPC contribution into ac¬ 
count in the px distributions. 

The Pz distributions are shown in Fig. Figure 
presents the ratios of LHCf Pz distributions to the Pz 
distributions predicted by the hadronic interaction mod¬ 
els. A similar tendency to that found in p-l-p collisions at 
y/s = 7 TeV is found for LHCf data relative to model pre¬ 
dictions. Concerning the comparison of hadronic interac¬ 
tion models with LHCf data, QGSJET shows a very good 
agreement at px < 0.2 GeV. However at px > 0.2 GeV, 
there are no models giving a consistent description of 
LHCf data within uncertainty over all Pz bins, although 
EPOS shows a certain compatibility with LHCf data for 
Px > 0.4 GeV and for pz < 3 TeV. The dpmjet pre¬ 
dictions agree with LHCf data at px < 0.6 GeV and 


Pz < 2 TeV, while showing a harder distribution at higher 
Pz similar to p -I- p collisions. Again note the character¬ 
istic bump found in the LHCf data at Pz 1.2 TeV and 
Px < 0.4 GeV, originating from the channel 7 +p —> 7r*^-|-p 
via baryon resonances in UPCs. 


VII. COMPARISONS OF THE LHCF 
MEASUREMENTS AMONG DIFFERENT 
COLLIDING HADRONS AND ENERGIES 

A. Average transverse momentum 

According to the scaling law proposed in Ref. jSS], the 
average transverse momentum, denoted (px)j as a func¬ 
tion of rapidity should be independent of the center-of- 
mass energy in the projectile fragmentation region. Here 
we obtain and compare the (px) distributions as func- 
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FIG. 9. (color online). LHCf Pt distributions (filled circles) in p + p collisions at = 2.76TeV. Error bars indicate 
the total statistical and systematic uncertainties. The predictions of hadronic interaction models are shown for comparison: 
DPMJET (solid red line), QGSJET (dashed blue line), sibyll (dotted green line), epos (dashed-dotted magenta line), and pythia 
(dashed-double-dotted brown line). 


tions of rapidity for p + p and p + Pb collisions. In the 
study of this paper, (px) is obtained by three methods 
discussed below. The first two methods use analytic dis¬ 
tributions that are fit to the LHCf data and the third 
uses numerical integration of the LHCf data. 

The first method uses the fit of an empirical Gaussian 
distribution to the LHCf px distributions for each rapid¬ 
ity range in Fig. mill and The second method uses 
a Hagedorn function. Here we pay attention to the fact 
that soft scattering dominates the measured 7 r° events for 
px <~ 1 GeV thus excluding from the analysis a power- 
law distribution that is used predominantly for hard scat¬ 
tering at px 1 GeV. These methods do not necessar¬ 
ily require that the measured px distribution be available 
down to 0.0 GeV, although the best-fit distribution may 
then include a systematic uncertainty in its fit | 66 | . De¬ 
tailed descriptions of the parametrization and derivation 
of (px) by using the best-fit Gaussian distribution can be 
found in Ref. m- In a Hagedorn function | 66 | . the in¬ 
variant cross section of identified hadrons, namely 7r°s in 
this paper, with a given mass m [GeV] and temperature 
T [GeV] can be written as 


CTinel dp-^ 


Pt 


fn^ Ki 

n—1 


,\/p| 


T 


( 6 ) 


where A [GeV“^[ is a normalization factor and Ki is the 
modified Bessel function. Approximately half of the 7 r° 
measured with the LHCf detector are daughters from 
the decay of parent baryons and mesons and are not di¬ 
rectly produced. Thus the measured px distribution is 
no longer a thermal distribution of prompt 7r°s and so 
we set m as a free parameter as well as A and T in the 
fit of a Hagedorn function to the px distribution. Equa¬ 
tion (|^ converges by n « 5 and the computation is in 
fact stopped at n = 10. The (px) value is calculated 
by using the modified Bessel functions AI 5/2 and K 2 as 
functions of the ratio of the best-fit m and T values [ 66 ] , 


(tt) = 


ITTmf K^/2{n{fh/T)) 
2 E“i^2(n(m/f)) ■ 


(7) 


For reference, Fig. 17 shows LHCf px distributions (filled 


black circles) and the best fits of the Gaussian distri¬ 
butions and the Hagedorn functions. The left panel of 
Fig. pT| shows the results for 9.2<?/<9.4inp-|-p colli¬ 
sions at ^/s = 7 TeV. The best-fit Gaussian distribution 
(dotted red curve) and Hagedorn function (dashed blue 
curve) to the LHCf data mostly overlap each other and 


give compatible (px) values. The right panel of Fig. [17 
shows the results for — 9.2 > piab > — 9.4 in p -|- Pb colli¬ 
sions at Gsnn = 5.02 TeV. Note that the px distribution 
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FIG. 10. (color online). Ratios of LHCf px distributions to the px distributions predicted by hadronic interaction models in 
p + p collisions at ^ = 2.76TeV are shown by solid red line (dpmjex), dashed blue line (qgsjex), dotted green line (sibyll), 
dashed-dotted magenta line (epos), and dashed-double-dotted brown line (pyxhia). Shaded areas indicate the range of total 
uncertainties of the px spectra. 
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FIG. 11. (color online). LHGf Pz distributions (filled circles) in p-f p collisions at ^/s = 2.76 TeV. Error bars indicate 
the total statistical and systematic uncertainties. The predictions of hadronic interaction models are shown for comparison: 
DPMJEX (solid red line), QGSJET (dashed blue line), sibyll (dotted green line), epos (dashed-dotted magenta line), and pyxhia 
(dashed-double-dotted brown line). 













































































17 




DPMJET 3.06 
EPOS LHC 
QGSJET 11-04 
SIBYLL 2.1 
PYTHIAS.185 


FIG. 12. (color online). Ratios of LHCf Pz distributions to the Pz distributions predicted by hadronic interaction models in 
p + p collisions at ^/s = 2.76TeV are shown by solid red line (dpmjet), dashed blue line (qgsjet), dotted green line (sibyll), 
dashed-dotted magenta line (epos), and dashed-double-dotted brown line (pythia). Shaded areas indicate the range of total 
uncertainties of the pz distributions. 


for LHCf data is plotted after subtraction of the UPC 
component where the systematic uncertainty in the sim¬ 
ulation of UPC events has been taken into account. The 
best-fit Gaussian distribution and the Hagedorn function 
reproduce the LHCf px distributions within the total un¬ 
certainties and are also compatible with each other. 

Finally, for the third method, (px) is obtained by nu¬ 
merically integrating the px distributions in Fig. min 
and The LHCf px distributions in p -h Pb colli¬ 
sions have already had the UPC component subtracted. 
In this approach, (px) is calculated only in the rapid¬ 
ity range where the px distributions are available down 
to 0.0 GeV. The high-px tail that extends beyond the 
data (px ^ (pt)) has a negligible contribution to (px). 
The final (px) values obtained in this analysis, denoted 
(PT)LHCf> have been determined by averaging the (px) 
values calculated with the three above described inde¬ 
pendent approaches: Gaussian, Hagedorn and numerical 
integration. The uncertainty of (px)LHCf for each ra¬ 
pidity bin is assigned to fully cover the minimum and 
maximum (px) values obtained by the three approaches. 
The (px)LHCf values are summarized in Table. |HI[ 

In Fig. 18 (px) in p -f p collisions at ^/s = 2.76 and 
Pb collisions at 


7TeV, and in p + Pb collisions at y^SNN = 5.02 TeV are 
presented as a function of rapidity loss Ay = pbeam — y, 
where pbeam is the beam rapidity for each collision en¬ 
ergy. The shift in rapidity by j/beam allows a direct com¬ 
parison to be made between the (px) results at different 
collision energies. We see that for Ay > —1.3, (px) at 
i/s = 2.76 TeV (open red circles) has slightly smaller val¬ 
ues than at y/s = 7TeV (filled black circles), although 
the two sets of data are mostly compatible at the ±10% 


level. For reference, the SppS UA7 results for p±p colli¬ 
sions at y/s = 630 GeV (open magenta squares) show 
a rapid roll off of (px) at low Ay compared to LHCf data. 
The LHCf and UA7 results are particularly incompatible 
for —0.3 < Ay < 0.3. The comparison of the LHCf 
data with the UA7 results indicates that (px) may de¬ 
pend on the center-of-mass energy. However, in order 
to firmly confirm a center-of-mass energy dependence of 
(px), we need to have experimental data at a lower col¬ 
lision energy, e.g., ^/s < ITeV and with a wider range 
of rapidity. Approved plans are underway to obtain this 
data by installing the LHCf detector at the RHIC ZDC 
position [^. The (px) values obtained fromp±Pb colli¬ 
sions at .^snn = 5.02 TeV (filled blue triangles) are con¬ 
sistent with those from p + p collisions at y/s = 7 TeV 
within the systematic uncertainties present. The pre¬ 
dictions from DPMJET (thick solid red line) and QGSJET 
(thin solid blue line) in p+p collisions at ^/s = 7 TeV and 
p±Pb collisions at ^^snn = 5.02 TeV have been added to 
for reference. The predictions at y/s = 2.76 TeV 
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Fig._ 

are excluded in Fig. since these curves mostly over¬ 
lap with those at 7 TeV. LHCf data in p ± p collisions 
at y/s = 7 TeV are close to the predictions made by dp¬ 
mjet at large Ay (small y) and become close to those 
made by QGSJET at small Ay (large y). These relations 
between LHCf data and the model predictions are con¬ 
sistent with the Px distributions shown in Fig. and 
The prediction from dpmjet (thick dashed red line) for 
p ± Pb collisions at .^snn = 5.02 TeV is compatible with 
the LHCf result for —0.3 < Ay < 0.2, which is derived 
from the good agreement of this model with LHCf data at 
—8.8 > piab > —10.0 and px < 0.3 GeV. Conversely, the 
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FIG. 13. (color online). LHCf pt distributions (filled circles) in p-|- Pb collisions at y^SNN = 5.02TeV. Error bars indicate 
the total statistical and systematic uncertainties. The predictions of hadronic interaction models are shown for comparison: 
DPMJET (solid red line), qgsjet (dashed blue line), and epos (dashed-dotted magenta line). 


prediction obtained from QGSJET (thin dashed blue line) 
is smaller than LHCf data for Ay > —0.5 and approaches 
the LHCf results on decreasing Ay. This tendency was 
already found in Fig. [T^ the prediction from QGSJET 
shows an overall agreement with LHCf pt distributions 
at yiab < -9.8. 


B. Limiting fragmentation 

The hypothesis of limiting fragmentation |121[I^ as¬ 
serts that the number of fragments of a colliding hadron 
will follow a limiting rapidity distribution in the rest 
frame of the target hadron. In this case the rapidity 
distribution of the secondary particles in the forward ra¬ 
pidity region would be independent of the center-of-mass 
energy. In this paper, a test of the limiting fragmenta¬ 


tion hypothesis is performed by using LHCf data in p-\-p 
collisions at y/s = 2.76 and 7TeV. 

The normalized rapidity distribution of 7r*^s, 
{1 /ainei) (da/dy), in this analysis can be obtained 
by using very similar methods that were used for the 
derivation of the average px in Sec. |VH A[ 

The first method uses the fit of an empirical distribu¬ 
tion to the LHCf Pt distributions in Fi g. [5] an d in each 
rapidity range. As we discussed in Sec. |VH^ two distri¬ 
butions are chosen to parametrize the px distributions: 
a Gaussian distribution and a Hagedorn function. The 
rapidity distribution is derived by integrating the best-fit 
Gaussian distribution and Hagedorn function along the 
Px axis from 0.0 GeV to infinity. 

The rapidity distribution can also be obtained by nu¬ 
merically integrating the px distributions in Fig.j^an m 
In this approach, the derivation of the (I/crinei)(do’/dy) 
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FIG. 14. (color online). Ratios of LHCf pt distributions to the pt distributions predicted by hadronic interaction models in 
p + Pb collisions at ^snn = 5.02 TeV are shown by solid red line (dpmjet), dashed blue line (qgsjet), and dashed-dotted 
magenta line (epos). Shaded areas indicate the range of total uncertainties of the px distributions. 


value is possible only in the rapidity range where the px 
distributions are available down to 0.0 GeV. Again, the 
final rapidity distribution is derived by averaging the ra¬ 
pidity distributions obtained by the above three methods. 
The estimated uncertainty is obtained from the minimum 
and maximum values for each rapidity bin. 


Figure shows the rapidity distributions as functions 
of the rapidity loss Ay (i.e., yheam — y) in P+P collisions at 
^/s = 2.76 TeV (open red circles) and 7 TeV (filled black 
circles). The rapidity distributions for both collision en¬ 
ergies mostly appear to lie along a common curve in the 
rapidity range —1.8 < Ay < —0.8. LHCf data are con¬ 
sistent at the ±15 % level with the hypothesis of limiting 
fragmentation in the very forward region. 

For comparison the experimental results from the UA7 
experiment m are also shown in Fig. |19[ The extrap¬ 


olated LHCf curve at 7TeV to higher Ay (i.e., lower y) 
could be compatible with the UA7 results, at least for 
Ay < 0.5. 


The predictions of dpmjet (thick red curve) and 
QGSJET (thin blue curve) at ^/s = 7 TeV have been 
added to Fig. 19 for reference. The predictions at y/s = 
2.76 TeV have been omitted, since these curves mostly 
overlap with those at 7 TeV since limiting fragmentation 
holds in dpmjet and QGSJET. The best agreement with 
LHCf data at ^/s = 2.76 and 7 TeV is obtained by the 
QGSJET model. The dpmjet predictions generally give 
a larger ir^ yield and a harder px distribution especially 
for y > 9.8 at = 7TeV and for y > 9.4 at 2.76 TeV. 
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FIG. 15. (color online). LHCf distributions (filled circles) in p + Pb collisions at ^snn = 5.02 TeV. Error bars indicate 
the total statistical and systematic uncertainties. The predictions of hadronic interaction models are shown for comparison: 
DPMJET (solid red line), qgsjet (dashed blue line), and epos (dashed-dotted magenta line). 


C. Feynman scaling 

In Ref. |17| Feynman proposed that the production 
cross sections of secondary particles as functions of the 
Feynman-x variable (defined by xp = ^') were in¬ 

dependent of the incident energy in the forward region. 
If the so-called Feynman scaling holds, the differential 
cross section as a function oi Xp (hereafter xp distribu¬ 
tion) (xpja\-ae\){d(Jjdxp) should be independent of the 
center-of-mass energy for xp > 0.2. Here the rapidity 
distribution introduced in Sec. IVII^ can be rewritten as 

\ da E da xp da 

^inel dy rTjnel dpz f^inel dxp 

where xp = 2E j^/s and dy = dpz/E are used for the 
second form. Considering p^ ^ E in the forward region, 
Xp can be considered as xp and thus the right hand side 
of Eq. (|^ becomes approximately {xp/ainei)ida/dxp). 
Consequently, the limiting fragmentation hypothesis that 
states {l/ainei)ida/dy) is independent of the center-of- 
mass energy in each rapidity bin can be rewritten as 


Feynman scaling which states {xp/ainei)(da/dxp) is in¬ 
dependent of the center-of-mass energy in each xp bin. 
In this paper, we test the Feynman scaling hypothesis by 
comparing LHCf data in p + p collisions at ^/s = 2.76 
and 7 TeV. 

In Fig. |20[ we compare the xp distributions in the px 
range 0.0 < px < 0.4 GeV. Other px ranges are excluded 
from the comparison, since LHCf data at ^/s = 2.76 TeV 
are unavailable outside this range. The xp distributions 
at y/s = 2.76 and 7TeV are compatible with each other 
at the ±20 % level. In Fig. we further compare the 
Xp distributions for the reduced px ranges: 0.0 < px < 
0.2 GeV and 0.2 < px < 0.4 GeV. At 0.0 < px < 0.2 GeV, 
only the bin 0.73 < xp < 0.82 at ^/s = 2.76 TeV deviates 
from the one at 7 TeV by 30%, while the other bins are 
consistent within their uncertainties. At 0.2 < px < 
0.4 GeV, all bins at ^/s = 2.76 TeV are consistent with 
the ones at 7 TeV, except for the bin 0.82 < xp < 0.91 
that has a smaller (40%) cross section than at 7 TeV, 
although there is a large uncertainty at 2.76 TeV. Overall 
the Xp distributions at ^/s = 2.76 and 7 TeV indicate 
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FIG. 16. (color online). Ratios of LHCf Pz distributions to the Pz distributions predicted by hadronic interaction models in 
p + Pb collisions at ^snn = 5.02 TeV are shown by solid red line (dpmjet), dashed blue line (qgsjet), and dashed-dotted 
magenta line (epos). Shaded areas indicate the range of total uncertainties of the px distributions. 


TABLE III. The average tt® transverse momenta for the ra¬ 
pidity range 8.8 < y < 10.6 in p -|-p collisions at -^i = 2.76 
and 7 TeV and for the rapidity range —8.8 > i/iab > —10.6 in 
p + Pb collisions at y'sNN = 5.02 TeV. 


Rapidity n (pT)LHCf [MeV] 

"p -i- p 2.76 TeV p-PpTTeV p -h Pb 5.02 TeV 


[8.8, 9.0] 

103.5 

± 

7.5 

242.8 

± 

8.6 

244.5 

± 

43.2 

[9.0, 9.2] 

78.5 

± 

7.8 

208.5 

± 

6.1 

223.1 

± 

12.7 

[9.2, 9.4] 

76.4 

± 

5.7 

182.6 

± 

4.3 

189.9 

± 

7.6 

[9.4, 9.6] 

60.3 

± 

5.2 

160.2 

± 

3.8 

173.8 

± 

17.2 

[9.6, 9.8] 

50.4 

± 

10.4 

132.3 

± 

3.4 

138.1 

± 

18.7 

[9.8, 10.0] 




113.9 

± 

3.4 

113.0 

± 

6.3 

[10.0,10.2] 




87.3 

± 

3.9 

112.2 

± 

15.4 

[10.2,10.4] 




67.5 

± 

3.0 

90.7 

± 

6.7 

[10.4,10.6] 




55.6 

± 

3.1 

61.0 

± 

6.6 


^ The rapidity values for p + Pb collisions are in the detector 
reference frame and must be multiplied by -1. 


that Feynman scaling holds at the ±20% level at these 
center-of-mass energies in the very forward region. 

Besides a test of the Feynman scaling, we find in Fig.[^ 
that the yield of 7r°s at i/s = 2.76 TeV relative to 7 TeV is 
slightly larger for 0.0 < pt < 0.2 GeV and slightly smaller 
for 0.2 < px < 0.4 GeV. This tendency means that the 
Pt distributions at y/s = 2.76 TeV are softer than those 
at 7TeV, leading to the small (px) values at 2.76 TeV 
relative to those at 7 TeV as already found in Fig. [T^ 


D. Pt dependence of the xp distributions 

In hadronic interactions at large rapidities, partons 
from the projectile and target hadrons generally have 
large and small momentum fractions respectively, since 
the momentum fraction that the parton itself carries rela¬ 
tive to the parent projectile and target hadrons, i.e., the 
Bjorken-a; variable or Xsj, is proportional to {+y 
for projectile and —y for target). Here we note that 
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FIG. 17. (color online). LHCf px distributions (filled black 
circles), the best-fit Gaussian distributions (dotted red curve), 
and the best-fit Hagedorn functions (dashed blue curve). Left: 
the data for p + p collisions at ^/s = 7TeV. Right: the data 
for p + Pb collisions at .^snn = 5.02 TeV. 
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FIG. 18. (color online). Average px as a function of rapidity 
loss Ap = pbeam — V- Open red circles and filled black cir¬ 
cles indicate LHGf data in p -f p collisions at ^/8 = 2.76 and 
7 TeV, respectively. The results of the UA7 experiment (open 
magenta box) at SppS (p + p collisions at ^/s = 630 GeV) and 
the predictions from dpmjex (thick lines) and qgsjex (thin 
lines) are added for reference. 


FIG. 19. (color online). The yield in each rapidity interval 
as a function of rapidity loss Ay = Pbeam —P- Open red circles 
and filled black circles indicate LHCf data in p -f p collisions 
at y's = 2.76 and 7 TeV, respectively. The results of the UA7 
experiment (open magenta squares) at SppS (p + p collisions 
at y'i = 630 GeV) and the predictions by dpmjex (thick red 
line) and qgsjex (thin blue line) are added for reference. 



a parton (dominantly gluon) density, rapidly increases 
with decreasing xbj when xbj < 0.01 with the target 
approaching the blackbody limit where the gluon den¬ 
sity is saturated. In the blackbody regime, the partons 
cannot go through the target nuclear medium without 
interaction and suffer transverse momenta transfers pro¬ 
portional to the saturation momentum scale Qs- The Qs 
values in the very forward region are ~ 1 GeV in p -|- p 
collisions and ~ 10 GeV in p-|-Pb collisions, although the 
calculation of Qs itself suffers from both theoretical and 
experimental uncertainties and is also dependent on the 
impact parameter of the colliding hadrons Ca¬ 
in the pt region below Qs, the xp distribution in the 


FIG. 20. (color online). The tt® yield at 0.0 < px < 0.4GeV 
as a function of x_f. Open red circles and filled black circles 
indicate LHCf data in p+p collisions at y/s = 2.76 and 7 TeV, 
respectively. 


forward region can be asymptotically written |69| as 


xp da 

^inel dXp 


OC (1 — Xp)°‘. 


(9) 


where a is the leading exponent. In the blackbody 
regime, the xp distribution of the leading hadron is 
strongly suppressed and thus a increases relative to the 
value found for a dilute target. Conversely, a decreases 
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FIG. 21. (color online). The 7r° yield in each px range as 
a function of xf- Left: the distributions for 0.0 < px < 
0.2 GeV. Right: the distributions for 0.2 < px < 0.4 GeV. 
Open red circles and filled black circles indicate LHGf data in 
p + p collisions at = 2.76 and 7TeV, respectively. 


with increasing px when px approaches or exceeds the 
saturation momentum scale Qs- 

Figurej^shows the best-fit leading exponent a in each 
Px range in p -|- p and p -|- Pb collisions. The leading ex¬ 
ponent in p -b p collisions at ^/s = 7 TeV (filled black 
circles) is a « 3.7 at px < 0.6 GeV and decreases to 
a « 3.0 at 0.6 < px < 1.0 GeV. The reduction of a with 
increasing px can be understood as much of the target 
staying in the blackbody regime for px < 0.6 GeV and 
then gradually escaping from the blackbody regime for 
Px > 0.6 GeV. The leading exponent in p-l-p collisions at 
■y/s = 2.76 TeV (open red circles) is slightly smaller than 
that at 7 TeV though with large uncertainty. The com¬ 
parison between -y/s = 2.76 TeV and 7 TeV may indicate 
that the upper px limit of the measurement at 2.76 TeV 
is near the saturation momentum at 2.76 TeV and that 
the suppression due to the gluon density is weaker than 
at 7TeV, although the calculated Qs at 2.76 TeV is only 
slightly different from the Qs at 7 TeV. The leading ex¬ 
ponents in p -b Pb collisions at = 5.02 TeV (filled 

blue triangles) are rather flat along the px axis, within 
uncertainties that are generally larger than those in p-bp 
collisions. This may indicate that the saturation momen¬ 
tum in p -b Pb collisions is well above the measured px 
range and also that the xp distributions in p -b Pb colli¬ 
sions are suppressed relative to those for p-bp collisions. 


E. Nuclear modification factor 


The effects of high gluon density in the target are in¬ 
ferred from the comparison of the leading exponent a 
between in p -b p and p -b Pb collisions (see the preceding 
Sec. VII Db Here we introduce the nuclear modification 
factor that quantifies the px spectra modification caused 
by nuclear effects in p -b Pb collisions with respect to 
p-bp collisions. The nuclear modification factor Rpph is 



FIG. 22. (color online). The best-fit leading exponent of (1- 
xf) as a function of px. Open red circles and filled black 
circles indicate LHGf data in p -b p collisions at ^/s = 2.76 
and 7 TeV, respectively. Filled blue triangles indicate LHGf 
data in p -b Pb collisions at y^SNN = 5.02 TeV. 


defined as 


-RpPb = 


^inel /dp^ 

(iVcon)CTPel’ -E^clVPP/dpa ’ 


( 10 ) 


where /dp^ and Ed^a'^'^/dp^ are the inclusive 

cross sections of 7r° production in p -b Pb collisions at 
.^snn = 5.02 TeV and in p-bp collisions at y/s = 5.02 TeV, 
respectively. These cross sections are obtained from 
Eq. (§ and Eq. 0: with the subtraction of the expected 
UPG contribution applied to the cross section for p -b Pb 
collisions. The uncertainty in the inelastic cross section 
crPjgJ’ is estimated to be ±5 % [T^ . The average number 
of binary nucleon-nucleon collisions in a p-bPb collision, 
(.^coii) = 6.9, is obtained from MG simulations using the 
Glauber model m- The uncertainty of crP,gj/(iVcou) is es¬ 
timated by varying the parameters in the calculation with 
the Glauber model and is of the order of ±3.5 % [19]. Fi¬ 
nally the quadratic sum of the uncertainties in crPjgJ’ and 
^inei/(^cou) is added to RpPh- 

Since there is no LHGf data for p ± p collisions at ex¬ 
actly y/s = 5.02 TeV, Ed^a'^'^/dp^ is derived by scaling 
the Px distributions taken in p ± p collisions to other 
collision energies. The derivation follows three steps. 
First, the (px) at y/s = 5.02 TeV is estimated by in¬ 
terpolating the measured (px) values at 7TeV. The un¬ 
certainty of the interpolated (px) values is estimated to 
be ±10% according to the differences between the mea¬ 
sured (px) values at y/s = 2.76 and 7TeV for —1.7 < 
Ay < —0.8 (see Fig. [^. Second, the absolute nor¬ 
malization of the Px distribution value in each rapidity 
range for y/s = 5.02TeV, i.e., {l/aiy^ei){da/dy), is deter¬ 
mined by interpolating the rapidity distribution at 7 TeV 
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(see Fig. 191. The uncertainty of the absolute normal¬ 
ization is estimated to be ±15% according to the dis¬ 
cussion in Sec. IVIIBI and is taken into account in the 
interpolated normalization. Finally, the px distributions 
at ^/s = 5.02 TeV are produced by assuming that the 
Pt distribution follows either a Gaussian distribution or 
a Hagedorn function and by using the (px) values ob¬ 
tained in the first step and the normalization obtained 
in the second step. The difference of the px distribution 
produced using a Gaussian distribution or a Hagedorn 
function gives the systematic uncertainty. Note that the 
rapidity shift —0.465 explained in Sec. |IIIB] is also taken 
into account for the px distribution in p ± p collisions at 
V^ = 5.02 TeV. 

Figure 23 shows the nuclear modification factors i?ppb 
obtained from LHGf data and the predictions from the 
hadronic interaction models, dpmjet (solid red curve), 
QGSJET (dashed blue curve), and EPOS (dotted magenta 
curve). LHGf data show a strong suppression with Rpph 
equal to 0.3 at j/iab — 8.8 down to < 0.1 at piab ~ —10.8, 
although a large uncertainty is present due to system¬ 
atic uncertainties in the estimation of the (px) values 
in p ± p collisions at ^/s = 5.02 TeV. All hadronic in¬ 
teraction models, which employ different approaches for 
the nuclear effects, predict small values of i?ppb < 0.15. 
Within the uncertainties the hadronic interaction mod¬ 
els show an overall good agreement with Rpph estimated 
from LHGf data. 


hold in the forward region at the 15%-20% level. The 
leading exponent a and the nuclear modification factor 
i?pPb derived from LHGf data indicate a strong suppres¬ 
sion of TT^ production from the nuclear target relative to 
that from the nucleon target. Within the uncertainties 
all of the hadronic interaction models presented gave an 
overall good agreement with Rppb estimated by LHGf 
data. According to the analysis in this paper, we expect 
that the number of particles leading to an electromag¬ 
netic component in air showers would follow the limit¬ 
ing rapidity distribution and Feynman scaling hypothe¬ 
ses. Gombining the results for forward 7 r°s in this paper 
with the recent results for forward neutrons in Ref. m 
strongly constrain models for air shower production at 
the TeV scale. 

As a future prospect, additional analyses using cor¬ 
relations between forward ir^s and other particles (e.g., 
two-particle angular correlations) are needed to reach a 
better understanding of the forward meson production 
mechanism and the strong suppression of production 
in p±Pb collisions compared to p+p collisions. The AT¬ 
LAS and LHGf Gollaborations have taken p + p data at 
= 13 TeV and p ± Pb data at y^SNN = 5.02 TeV with 
common triggers. This data could provide the possibility 
for performing analyses of two particle correlations. 
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APPENDIX: DATA TABLES 


The inclusive production rates of 7 r°s measured by 
LH Gf w i th all c orrections applied are summarized in Ta¬ 
bles IV - XVHI The ratios of 7 r° production rate of MG 
simulation to data are summarized in Tables IXIXLILVI 
The LHGf px distributions for p ± Pb collisions have the 
UPG component subtracted. The nuclear modification 
factor of 7 r°s obtained from LHGf data are summarized 


in Tables iLVD-rnTXl 
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Pt [GeV] [GeV] 

FIG. 23. (color online). The nuclear modification factor for tt^s. Filled circles indicate the factors obtained from LHCf data. 
Error bars indicate the total uncertainties incorporating both statistical and systematic uncertainties. Other lines are the 
predictions from hadronic interaction models: dpmjet (solid red line), qgsjet (dashed blue line), and epos (dashed-dotted 
magenta line). 
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TABLE IV. Production rate for the tt® production in the rapidity range 8.8 < y < 9.0 in p + p collisions and in the rapidity 
range —8.8 > yiab > —9.0 in p + Pb collisions. The rate and corresponding total uncertainty are in units of [GeV“^]. 


Pt range [GeV] 

p + p at i/i = 

2.76 TeV 

p + p at 

= 7 TeV 

p ± Pb at ,/snn 

= 5.02 TeV 

Production rate 

Total error 

Production rate 

Total error 

Production rate 

Total error 

[0.10,0.15] 

3.71x10''^ 

±8.99x10“^ 

2.56x10“^ 

±3.75x10“® 

1.63x10“® 

±5.65x10“® 

[0.15,0.20] 

1.92x10"^ 

±3.36x10"^ 

1.98x10“^ 

±2.05x10“® 

1.86x10“® 

±4.19x10“® 

[0.20,0.25] 

3.47x10"^ 

±7.22x10-^ 

1.20x10“^ 

±1.29x10“® 

1.52x10“® 

±2.76x10“® 

[0.25,0.30] 

5.67x10"^ 

±2.49x10"^ 

6.96x10“^ 

±5.65x10“® 

1.06x10“® 

±1.52x10“® 

[0.30,0.35] 



3.94x10“^ 

±3.43x10“® 

5.59x10“® 

±8.43x10“® 

[0.35,0.40] 



2.23x10“^ 

±2.28x10“® 

3.80x10“® 

±6.14x10“® 

[0.40,0.45] 



1.26x10“^ 

±1.52x10“® 

2.13x10“® 

±3.86x10“® 

[0.45,0.50] 



7.06x10“® 

±8.82xl0“'‘ 

9.51x10“® 

±1.84x10“® 

[0.50,0.55] 



4.21x10“® 

±4.93xl0“'‘ 

6.56x10“® 

±1.09x10“® 

[0.55,0.60] 



2.31x10“® 

±2.61xl0“'‘ 

3.18x10“® 

±5.64x10“® 

[0.60,0.65] 



1.15x10“® 

±1.44xl0“'‘ 

1.64x10“® 

±3.28x10“® 

[0.65,0.70] 



4.36xl0“'‘ 

±7.06x10“® 

8.98x10“"® 

±1.95x10“® 

[0.70,0.75] 



2.12xl0“'‘ 

±4.56x10“® 

5.88x10“'® 

±1.60x10“® 

[0.75,0.80] 



1.38xl0“'‘ 

±3.15x10“® 

1.54x10“® 

±7.96x10“® 

[0.80,0.85] 



3.98x10“® 

±1.59x10“® 

3.63x10“® 

±4.26x10“® 

[0.85,0.90] 



4.52x10“® 

±3.26x10“® 

7.07x10“® 

±4.06x10“® 

[0.90,0.95] 



7.09x10“® 

±4.54x10“® 

< 1.48x10“® 



TABLE V. Production rate for the tt® production in the rapidity range 9.0 < y < 9.2 in p + p collisions and in the rapidity 
range —9.0 > piab > —9.2 in p + Pb collisions. The rate and corresponding total uncertainty are in units of [GeV“^]. 


Pt range ]GeV] 

p ± p at i/s = 

2.76 TeV 

p ± p at i/s 

= 7 TeV 

p ± Pb at ,/snn 

= 5.02 TeV 

Production rate 

Total error 

Production rate 

Total error 

Production rate 

Total error 

[0.05,0.10] 

4.26x10“® 

±1.08x10“® 

2.98x10“® 

±4.72x10“® 

1.66x10“® 

±5.08x10“® 

[0.10,0.15] 

2.18x10“® 

±3.74x10“® 

2.27x10“® 

±2.58x10“® 

1.78x10“® 

±3.66x10“® 

[0.15,0.20] 

4.20x10“® 

±7.14x10“® 

1.42x10“® 

±1.12x10“® 

1.69x10“® 

±2.97x10“® 

[0.20,0.25] 

2.89x10“® 

±1.59x10“® 

8.48x10“® 

±6.20x10“® 

1.14x10“® 

±2.10x10“® 

[0.25,0.30] 



5.50x10“® 

±4.43x10“® 

5.31x10“® 

±8.70x10“® 

[0.30,0.35] 



3.28x10“® 

±3.02x10“® 

3.46x10“® 

±5.49x10“® 

[0.35,0.40] 



1.53x10“® 

±1.59x10“® 

1.73x10“® 

±3.28x10“® 

[0.40,0.45] 



6.83x10“® 

±8.38x10“® 

9.01x10“® 

±2.06x10“® 

[0.45,0.50] 



3.31x10“® 

±4.19x10“® 

4.87x10“® 

±1.15x10“® 

[0.50,0.55] 



1.35x10“® 

±1.74x10“® 

2.96x10“® 

±6.76x10“® 

[0.55,0.60] 



7.21x10“® 

±8.11x10“® 

6.98x10“® 

±2.73x10“® 

[0.60,0.65] 



1.94x10“® 

±4.84x10“® 

3.24x10“® 

±1.66x10“® 

[0.65,0.70] 



9.56x10“® 

±2.79x10“® 

2.92x10“® 

±1.13x10“® 

[0.70,0.75] 



1.26x10“® 

±5.76x10“® 

1.56x10“® 

±4.63x10“® 

[0.75,0.80] 



5.05x10“® 

±3.11x10“® 

1.60x10“® 

±2.68x10“® 



28 


TABLE VI. Production rate for the tt® production in the rapidity range 9.2 < y < 9.4 in p + p collisions and in the rapidity 
range —9.2 > yiab > —9.4 in p + Pb collisions. The rate and corresponding total uncertainty are in units of [GeV“^]. 


Pt range [GeV] 

p + p at = 

2.76 TeV 

p ± p at 

= 7 TeV 

p ± Pb at ,/snn 

= 5.02 TeV 

Production rate 

Total error 

Production rate 

Total error 

Production rate 

Total error 

[0.00,0.05] 

4.46 xlO""^ 

±1.33x10“'' 

2.74x10“^ 

±5.08x10“® 

1.46x10“® 

±5.35x10“® 

[0.05,0.10] 

2.07x10"^ 

±4.12x10“® 

1.96x10“® 

±3.22x10“® 

1.71x10“® 

±3.98x10“® 

[0.10,0.15] 

7.15x10"^ 

±1.12x10“® 

1.54x10“® 

±1.76x10“® 

1.52x10“® 

±2.67x10“® 

[0.15,0.20] 

8.30x10"^ 

±2.24x10“"^ 

1.01x10“® 

±7.24x10“® 

1.39x10“® 

±2.42x10“® 

[0.20,0.25] 

4.11x10“® 

±4.14x10“® 

5.53x10“® 

±4.20x10“® 

6.48x10“® 

±1.15x10“® 

[0.25,0.30] 



3.26x10“® 

±2.77x10“® 

2.80x10“® 

±4.62x10“® 

[0.30,0.35] 



1.59x10“® 

±1.53x10“® 

1.57x10“® 

±2.87x10“® 

[0.35,0.40] 



7.26x10“® 

±9.04x10“"® 

8.06x10“® 

±2.12x10“® 

[0.40,0.45] 



1.72x10“® 

±3.83x10“"® 

4.43x10“® 

±1.46x10“® 

[0.45,0.50] 



4.12x10“"® 

±1.78x10“"® 

1.94x10“® 

±9.94x10“® 

[0.50,0.55] 



9.38x10“® 

±9.38x10“® 

1.48x10“"® 

±3.75x10“® 

[0.55,0.60] 



9.53x10“® 

±8.49x10“® 

1.28x10“® 

±2.08x10“® 


TABLE VII. Production rate for the rr® production in the rapidity range 9.4 < y < 9.6 in p + p collisions and in the rapidity 
range —9.4 > yiab > —9.6 in p + Pb collisions. The rate and corresponding total uncertainty are in units of [GeV“^]. 


Pt range ]GeV] 

p ± p at i/s = 

2.76 TeV 

p ± p at i/s 

= 7 TeV 

p ± Pb at ,/snn 

= 5.02 TeV 

Production rate 

Total error 

Production rate 

Total error 

Production rate 

Total error 

[0.00,0.05] 

1.71x10“® 

±4.08x10“® 

1.95x10“® 

±3.27x10“® 

1.39x10“® 

±3.89x10“® 

[0.05,0.10] 

6.51x10“® 

±1.29x10“® 

1.52x10“® 

±2.23x10“® 

1.36x10“® 

±2.90x10“® 

[0.10,0.15] 

7.27x10“® 

±1.97x10“® 

1.08x10“® 

±8.30x10“® 

1.17x10“® 

±2.22x10“® 

[0.15,0.20] 

3.84x10“® 

±3.87x10“® 

6.08x10“® 

±4.41x10“® 

7.42x10“® 

±1.35x10“® 

[0.20,0.25] 



3.37x10“® 

±2.85x10“® 

3.08x10“® 

±4.97x10“® 

[0.25,0.30] 



1.49x10“® 

±1.49x10“® 

1.55x10“® 

±2.91x10“® 

[0.30,0.35] 



5.36x10“® 

±7.82x10“® 

1.20x10“® 

±2.53x10“® 

[0.35,0.40] 



1.60x10“® 

±4.36x10“® 

4.64x10“® 

±1.83x10“® 

[0.40,0.45] 



2.72x10“® 

±1.45x10“® 

9.01x10“® 

±7.91x10“® 


TABLE VIII. Production rate for the tt® production in the rapidity range 9.6 < y < 9.8 in p + p collisions and in the rapidity 
range —9.6 > yiab > —9.8 in p + Pb collisions. The rate and corresponding total uncertainty are in units of [GeV“^]. 


Pt range ]GeV] 

p ± p at \/i = 

2.76 TeV 

p ± p at i/i 

= 7 TeV 

p ± Pb at y^iNN 

= 5.02 TeV 

Production rate 

Total error 

Production rate 

Total error 

Production rate 

Total error 

[0.00,0.05] 

3.56x10“® 

±1.06x10“® 

1.63x10“® 

±1.98x10“® 

1.17x10“® 

±3.27x10“® 

[0.05,0.10] 

7.59x10“® 

±2.56x10“® 

1.07x10“® 

±1.07x10“® 

1.09x10“® 

±2.36x10“® 

[0.10,0.15] 

4.71x10“® 

±4.75x10“® 

6.40x10“® 

±4.56x10“® 

6.56x10“® 

±1.41x10“® 

[0.15,0.20] 



3.51x10“® 

±2.77x10“® 

3.03x10“® 

±5.17x10“® 

[0.20,0.25] 



1.54x10“® 

±1.65x10“® 

1.42x10“® 

±2.67x10“® 

[0.25,0.30] 



3.27x10“® 

±6.32x10“® 

6.67x10“® 

±1.82x10“® 

[0.30,0.35] 



9.33x10“® 

±3.58x10“® 

1.22x10“® 

±8.18x10“® 

[0.35,0.40] 



1.69x10“® 

±1.30x10“® 

< 1.73x10“® 



TABLE IX. Production rate for the rr® production in the rapidity range 9.8 < y < 10.0 in p + p collisions and in the rapidity 
range —9.8 > yiab > —10.0 in p + Pb collisions. The rate and corresponding total uncertainty are in units of [GeV“^]. 


Pt range ]GeV] 

p ± p at y"s = 

2.76 TeV 

p ± p at y"s 

= 7 TeV 

p ± Pb at y'sNN 

= 5.02 TeV 

Production rate 

Total error 

Production rate 

Total error 

Production rate 

Total error 

[0.00,0.05] 



8.80x10“® 

±1.05x10“® 

9.81x10“® 

±3.22x10“® 

[0.05,0.10] 



5.90x10“® 

±5.65x10“® 

6.41x10“® 

±1.61x10“® 

[0.10,0.15] 



3.45x10“® 

±2.73x10“® 

2.42x10“® 

±5.51x10“® 

[0.15,0.20] 



1.45x10“® 

±1.37x10“® 

1.29x10“® 

±2.70x10“® 

[0.20,0.25] 



3.07x10“® 

±5.82x10“® 

7.63x10“® 

±1.84x10“® 

[0.25,0.30] 



5.83x10“® 

±2.37x10“® 

1.01x10“® 

±7.40x10“® 
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TABLE X. Production rate for the production in the rapidity range 10.0 < y < 10.2 in p + p collisions and in the rapidity 
range —10.0 > yub > —10.2 in p + Pb collisions. The rate and corresponding total uncertainty are in units of [GeV“^]. 


Pt range [GeV] 

p + p at = 

2.76 TeV 

p + p at i/i 

= 7 TeV 

p + Pb at ysNN = 5.02 TeV 

Production rate 

Total error 

Production rate 

Total error 

Production rate Total error 

[0.00,0.05] 



4.54xl0~^ 

±5.63x10"® 

4.38x10"® ±1.80x10"® 

[0.05,0.10] 



3.06x10"^ 

±2.95x10"® 

2.51x10"® ±6.11x10"® 

[0.10,0.15] 



1.24x10"^ 

±1.12x10"® 

7.25x10"® ±2.38x10"® 

[0.15,0.20] 



2.04x10"® 

±3.16x10""* 

7.17x10"® ±1.82x10"® 

[0.20,0.25] 



5.33x10"^ 

±2.41x10""* 

3.00x10""* ±5.46x10""* 


TABLE XL Production rate for the tt® production in the rapidity range 10.2 < y < 10.4 in p + p collisions and in the rapidity 
range —10.2 > i/iab > —10.4 in p + Pb collisions. The rate and corresponding total uncertainty are in units of [GeV“^]. 


Pt range ]GeV] 

p + p at = 

2.76 TeV 

p + p at 

= 7 TeV 

p + Pb at y'sNN = 5.02 TeV 

Production rate 

Total error 

Production rate 

Total error 

Production rate Total error 

[0.00,0.05] 



2.08x10"® 

±2.31x10"® 

6.38x10"® ±6.42x10"® 

[0.05,0.10] 



1.10x10"® 

±1.00x10"® 

9.20x10"® ±1.91x10"® 

[0.10,0.15] 



2.22x10"® 

±2.16x10"* 

5.33x10"® ±8.77x10"* 

[0.15,0.20] 



8.39x10"® 

±2.03x10"® 

4.87x10"* ±1.64x10"* 


TABLE XII. Production rate for the tt® production in the rapidity range 10.4 < y < 10.6 in p + p collisions and in the rapidity 
range —10.4 > piab > —10.6 in p + Pb collisions. The rate and corresponding total uncertainty are In units of [GeV“^]. 


Pt range ]GeV] 

p + p at = 

2.76 TeV 

p + p at 

= 7 TeV 

p + Pb at y'sNN = 5.02 TeV 

Production rate 

Total error 

Production rate 

Total error 

Production rate Total error 

[0.00,0.05] 



5.14x10"® 

±7.21x10"* 

8.01x10"® ±4.52x10"® 

[0.05,0.10] 



1.87x10"® 

±2.11x10"* 

1.78x10"® ±1.39x10"® 

[0.10,0.15] 



1.04x10"* 

±2.51x10"® 

1.04x10"* ±3.45x10"* 

[0.15,0.20] 



2.09x10"® 

±1.37x10"® 

2.82x10"* ±1.51x10"* 


TABLE XIII. Production rate for the tt® production in the rapidity range 10.6 < y < 10.8 in p + p collisions and in the rapidity 
range —10.6 > yiab > —10.8 in p + Pb collisions. The rate and corresponding total uncertainty are in units of [GeV“^]. 


Pt range ]GeV] 

p + p at = 

2.76 TeV 

p + p at i/s 

= 7 TeV 

p + Pb at y'sNN = 5.02 TeV 

Production rate 

Total error 

Production rate 

Total error 

Production rate Total error 

[0.00,0.05] 



6.81x10"* 

±1.68x10"* 

1.17x10"® ±1.18x10"® 

[0.05,0.10] 



1.42x10"* 

±3.84x10"® 

< 9.69x10"® 

[0.10,0.15] 



2.81x10"® 

±1.57x10"® 

2.85x10"® ±2.57x10"® 
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TABLE XIV. Production rate for the tt® production in the px range 0.0 < px < 0.2 GeV in p + p and p + Pb collisions. The 
rate and corresponding total uncertainty are in units of [GeV“^]. 


Pz range [GeV] 

p + p at x/i = 

2.76 TeV 

p ± p at i/i 

= 7 TeV 

p ± Pb at ^snn 

= 5.02 TeV 

Production rate 

Total error 

Production rate 

Total error 

Production rate 

Total error 

[ 625, 750] 

4.55x10"'' 

±1.85x10“" 





[ 750, 875] 

2.59x10“" 

±5.79x10“" 





[ 875,1000] 

1.19x10“" 

±1.85x10“" 

1.82x10“" 

±2.21x10“" 

1.59x10“" 

±3.25x10“" 

[1000,1125] 

3.83x10“'^ 

±5.87x10“^ 

1.52x10“" 

±1.47x10“" 

1.48x10“" 

±2.50x10“" 

[1125,1250] 

9.35xl0“'‘ 

±2.19x10“^ 

1.39x10“" 

±9.38x10“" 

1.48x10“" 

±2.41x10“" 

[1250,1380] 

1.16xl0“'‘ 

±7.08x10“® 

9.97x10“" 

±6.89x10“" 

1.21x10“" 

±2.18x10“" 

[1380,1500] 



7.54x10“" 

±5.69x10“" 

8.68x10“" 

±1.81x10“" 

[1500,1625] 



5.85x10“" 

±4.70x10“" 

5.79x10“" 

±1.31x10“" 

[1625,1750] 



4.71x10“" 

±3.99x10“" 

3.26x10“" 

±8.00x10“" 

[1750,1875] 



3.73x10“" 

±3.07x10“" 

2.42x10“" 

±5.41x10“" 

[1875,2000] 



2.84x10“" 

±2.45x10“" 

1.24x10“" 

±3.77x10“" 

[2000,2125] 



2.10x10“" 

±1.91x10“" 

1.60x10“" 

±3.58x10“" 

[2125,2250] 



1.42x10“" 

±1.34x10“" 

6.87x10“" 

±2.68x10“" 

[2250,2375] 



1.13x10“" 

±1.14x10“" 

7.45x10“" 

±2.37x10“" 

[2375,2500] 



5.48x10“" 

±6.88x10“" 

9.05x10“" 

±2.35x10“" 

[2500,2625] 



3.59x10“" 

±4.65x10“" 

6.81x10“" 

±1.77x10“" 

[2625,2750] 



2.38x10“" 

±3.28x10“" 

4.40x10“" 

±1.32x10“" 

[2750,2875] 



1.48x10“" 

±2.15x10“" 

3.88x10“" 

±9.52x10“^ 

[2875,3000] 



8.28x10“" 

±1.28x10“" 

3.44x10“" 

±7.94x10“^ 

[3000,3125] 



6.06x10“" 

±9.65x10“® 

1.83x10“" 

±5.02x10“^ 

[3125,3250] 



2.19x10“" 

±4.44x10“® 

< 9.58x10“® 


[3250,3375] 



4.76x10“® 

±2.06x10“® 

< 9.46x10“® 


[3375,3500] 



5.33x10“® 

±2.16x10“® 

< 2.35x10“® 


[3500,3625] 





< 3.85x10“® 


[3625,3750] 





< 5.17x10“® 


[3750,3875] 





< 1.46x10“® 


[3875,4000] 





< 1.18x10“® 
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TABLE XV. Production rate for the tt® production in the px range 0.2 < px < 0.4 GeV in p + p and p + Pb collisions. The 
rate and corresponding total uncertainty are in units of [GeV“^]. 


Pz range [GeV] 

p + p at i/s = 

2.76 TeV 

p ± p at i/s 

= 7 TeV 

p ± Pb at i/snn 

= 5.02 TeV 

Production rate 

Total error 

Production rate 

Total error 

Production rate 

Total error 

[ 750, 875] 

1.16x10"'' 

±3.73x10"® 





[ 875,1000] 

4.09x10“'^ 

±9.00x10"^ 



1.26x10"® 

±2.61x10"" 

[1000,1125] 

1.21x10“'^ 

±2.89x10"^ 

8.47x10"" 

±8.41x10"® 

1.25x10"® 

±1.78x10"" 

[1125,1250] 

5.55x10"® 

±5.59x10"® 

7.00x10"" 

±7.22x10"® 

9.80x10"" 

±1.44x10"" 

[1250,1380] 

3.38x10"® 

±3.41x10"® 

5.58x10"" 

±4.31x10"® 

7.55x10"" 

±1.15x10"" 

[1380,1500] 



4.55x10"" 

±3.76x10"® 

4.95x10"" 

±8.07x10"® 

[1500,1625] 



3.57x10"" 

±3.08x10"® 

3.39x10"" 

±5.72x10"® 

[1625,1750] 



3.05x10"" 

±2.78x10"® 

3.17x10"" 

±5.33x10"® 

[1750,1875] 



2.38x10"" 

±2.25x10"® 

2.35x10"" 

±4.11x10"® 

[1875,2000] 



1.80x10"" 

±1.88x10"® 

1.47x10"" 

±3.08x10"® 

[2000,2125] 



1.38x10"" 

±1.51x10"® 

1.15x10"" 

±2.64x10"® 

[2125,2250] 



8.76x10"® 

±1.14x10"® 

9.83x10"® 

±2.44x10"® 

[2250,2375] 



5.69x10"® 

±8.68x10"^ 

8.19x10"® 

±2.22x10"® 

[2375,2500] 



3.91x10"® 

±6.97x10"^ 

7.74x10"® 

±2.23x10"® 

[2500,2625] 



1.61x10"® 

±4.20x10"^ 

7.17x10"® 

±2.04x10"® 

[2625,2750] 



1.13x10"® 

±3.52x10"^ 

3.08x10"® 

±1.30x10"® 

[2750,2875] 



6.47x10-^ 

±2.13x10"^ 

2.05x10"® 

±9.47x10"® 

[2875,3000] 



4.19xl0"'‘ 

±1.42x10"^ 

1.95x10"® 

±7.89x10"® 

[3000,3125] 



2.21x10-^ 

±9.01x10"® 

8.21x10-® 

±5.73x10"® 

[3125,3250] 



8.45x10"® 

±4.42x10"® 

2.67x10"® 

±4.03x10"® 

[3250,3375] 



3.90x10"® 

±1.75x10"® 

3.56x10"® 

±3.27x10"® 

[3375,3500] 



1.75x10"® 

±1.14x10"® 

5.10x10"® 

±2.14x10"® 

[3500,3625] 





2.34x10"® 

±2.26x10"® 

[3625,3750] 





1.85x10"® 

±1.69x10"® 

[3750,3875] 





1.66x10"® 

±1.45x10"® 

[3875,4000] 





1.16x10"® 

±8.13x10"® 


TABLE XVI. Production rate for the tt® production in the px range 0.4 < px < 0.6 GeV in p + p and p + Pb collisions. The 
rate and corresponding total uncertainty are in units of [GeV“^]. 


Pz range [GeV] 

p ± p at i/i = 

2.76 TeV 

p ± p at i/i 

= 7 TeV 

p ± Pb at 

= 5.02 TeV 

Production rate 

Total error 

Production rate 

Total error 

Production rate 

Total error 

[1380,1500] 





2.68x10"" 

±5.38x10"® 

[1500,1625] 



1.15x10"" 

±3.18x10"® 

1.45x10"" 

±3.06x10"® 

[1625,1750] 



1.09x10"" 

±2.06x10"® 

1.49x10"" 

±2.48x10"® 

[1750,1875] 



9.86x10"® 

±9.93x10"® 

1.16x10"" 

±1.82x10"® 

[1875,2000] 



6.78x10"® 

±6.85x10"® 

8.25x10"® 

±1.27x10"® 

[2000,2125] 



4.66x10"® 

±3.96x10"® 

6.36x10"® 

±1.05x10"® 

[2125,2250] 



4.08x10"® 

±3.58x10"® 

4.98x10"® 

±9.15x10"® 

[2250,2375] 



2.44x10"® 

±2.38x10"® 

4.15x10"® 

±8.31x10"® 

[2375,2500] 



1.73x10"® 

±1.89x10"® 

1.69x10"® 

±6.44x10"® 

[2500,2625] 



1.12x10"® 

±1.86x10"® 

2.43x10"® 

±7.14x10"® 

[2625,2750] 



5.57x10"® 

±1.47x10"® 

1.24x10"® 

±6.12x10"® 

[2750,2875] 



2.93x10"® 

±1.07x10"® 

5.21x10"® 

±4.99x10"® 

[2875,3000] 



2.36x10"® 

±9.47x10"® 

1.37x10"® 

±6.40x10"® 
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TABLE XVII. Production rate for the tt® production in the px range 0.6 < px < 0.8 GeV in p + p and p + Pb collisions. The 
rate and corresponding total uncertainty are in units of [GeV“^]. 


Pz range [GeV] 

p + p at i/i = 

2.76 TeV 

p + p at i/i 

= 7 TeV 

p ± Pb at 

= 5.02 TeV 

Production rate 

Total error 

Production rate 

Total error 

Production rate 

Total error 

[1750,1875] 





4.14x10"® 

±1.08x10"® 

[1875,2000] 





3.80x10"® 

±6.61x10"® 

[2000,2125] 



2.47x10"® 

±3.18xl0"'‘ 

3.81x10"® 

±5.91x10"® 

[2125,2250] 



1.56x10"® 

±2.03xl0"'‘ 

2.24x10"® 

±3.91x10"® 

[2250,2375] 



1.26x10"® 

±1.57xl0"'‘ 

1.89x10"® 

±3.19x10"® 

[2375,2500] 



7.10xl0"'‘ 

±6.10x10"® 

1.47x10"® 

±2.64x10"® 

[2500,2625] 



5.30x10"^ 

±4.94x10"® 

9.45x10"® 

±2.05x10"® 

[2625,2750] 



3.82xl0"'‘ 

±4.62x10"® 

4.17x10"® 

±1.47x10"® 

[2750,2875] 



2.55xl0"'‘ 

±4.86x10"® 

4.56x10"® 

±1.35x10"® 

[2875,3000] 



1.20xl0"'‘ 

±3.30x10"® 

8.61x10"® 

±9.66x10"® 

[3000,3125] 



5.73x10"® 

±2.27x10"® 

< 3.48x10"® 


[3125,3250] 



3.21x10"® 

±1.74x10"® 

2.74x10"® 

±1.24x10"® 

[3250,3375] 



3.16x10"® 

±1.79x10"® 

1.15x10"® 

±8.87x10"® 

[3375,3500] 





2.97x10"® 

±5.38x10"® 


TABLE XVIII. Production rate for the tt® production in the px range 0.8 < px < 1.0 GeV 
rate and corresponding total uncertainty are in units of [GeV"®[. 

in p ± p and p ± Pb collisions. The 

Pz range [GeV] 

p ± p at = 2.76 TeV 

p ± p at 

= 7 TeV 

p ± Pb at ^snn 

= 5.02 TeV 

Production rate Total error 

Production rate 

Total error 

Production rate 

Total error 

[2250,2375] 


5.33x10"® 

±3.02x10"® 



[2375,2500] 


3.19x10"® 

±1.41x10"® 

5.46x10"® 

±2.57x10"® 

[2500,2625] 


2.08x10"® 

±7.78x10"® 

4.85x10"® 

±1.37x10"® 

[2625,2750] 


1.19x10"® 

±4.19x10"® 

2.69x10"® 

±1.04x10"® 

[2750,2875] 


1.11x10"® 

±3.70x10"® 

2.70x10"® 

±9.97x10"® 

[2875,3000] 


3.35x10"® 

±1.59x10"® 

1.84x10"® 

±7.34x10"® 

[3000,3125] 


2.98x10"® 

±1.29x10"® 

6.77x10"® 

±4.60x10"® 

[3125,3250] 


5.22x10"® 

±3.45x10"® 

< 1.10x10"® 


[3250,3375] 


3.32x10"® 

±2.53x10"® 



[3375,3500] 




< 1.23x10"® 
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TABLE XIX. Ratio of tt® production rate of MC simulation 
to data in the rapidity range 8.8 < y < 9.0 in p + p collisions 
at = 7 TeV. 

TABLE XXL Ratio of tt® production rate of MC simulation 
to data in the rapidity range 9.2 < y < 9.4 in p + p collisions 
at Vi = 7 TeV. 

pT range 
[GeV] 

DPMJET 

3.06 

QGSJET 

11-04 

SIBYLL 

2.1 

EPOS 

LHC 

PYTHIA 

8.185 

pT range dpmjet 

[GeV] 3.06 

QGSJET 

11-04 

SIBYLL 

2.1 

EPOS 

LHC 

PYTHIA 

8.185 

[0.10, 0.15] 

2.38 

2.09 

1.26 

1.89 

2.34 

[0.00, 0.05] 

1.53 

1.71 

0.72 

1.20 

1.46 

[0.15, 0.20] 

1.83 

1.43 

0.93 

1.43 

1.79 

[0.05, 0.10] 

1.41 

1.40 

0.64 

1.12 

1.34 

[0.20, 0.25] 

1.43 

0.97 

0.71 

1.02 

1.40 

[0.10, 0.15] 

1.23 

1.04 

0.52 

0.93 

1.15 

[0.25, 0.30] 

1.64 

0.94 

0.80 

1.15 

1.60 

[0.15, 0.20] 

1.18 

0.84 

0.50 

0.76 

1.12 

[0.30, 0.35] 

1.94 

0.97 

0.96 

1.40 

1.85 

[0.20, 0.25] 

1.39 

0.78 

0.56 

0.95 

1.24 

[0.35, 0.40] 

2.31 

1.02 

1.19 

1.45 

2.15 

[0.25, 0.30] 

1.53 

0.72 

0.60 

1.00 

1.28 

[0.40, 0.45] 

2.81 

1.03 

1.49 

1.43 

2.52 

[0.30, 0.35] 

2.06 

0.71 

0.76 

0.88 

1.59 

[0.45, 0.50] 

3.44 

0.94 

1.83 

1.47 

2.93 

[0.35, 0.40] 

2.91 

0.56 

0.94 

0.88 

2.07 

[0.50, 0.55] 

3.91 

0.78 

2.08 

1.48 

3.19 

[0.40, 0.45] 

7.74 

0.80 

2.17 

1.95 

4.87 

[0.55, 0.60] 

4.72 

0.67 

2.52 

1.66 

3.70 

[0.45, 0.50] 

18.99 

1.13 

4.27 

3.54 

10.20 

[0.60, 0.65] 

6.24 

0.65 

3.28 

1.88 

4.60 

[0.50, 0.55] 

40.53 

1.21 

7.25 

7.54 

17.05 

[0.65, 0.70] 

10.40 

0.81 

5.33 

2.77 

7.30 

[0.55, 0.60] 

13.68 

0.32 

2.01 

3.24 

4.50 

[0.70, 0.75] 

12.58 

0.73 

6.29 

3.14 

8.37 







[0.75, 0.80] 

10.66 

0.46 

4.96 

2.95 

6.69 







[0.80, 0.85] 

17.83 

0.45 

7.65 

5.60 

11.19 







[0.85, 0.90] 

62.69 

1.50 

21.76 

27.76 

35.11 







[0.90, 0.95] 

11.44 

0.24 

3.64 

8.26 

7.36 







[0.95, 1.00] 

16.38 

0.20 

4.59 

26.72 

16.18 













TABLE XXII. 

Ratio of 

TT® production rate of MC simulation 







to data in the rapidity range 9.4 < y < 9.6 inp + p collisions 







at Vi = 7 TeV. 






TABLE XX. Ratio of tt 

^ production rate of MC simulation to 







data in the rapidity range 9.0 < y < 9.2 in p + p collisions at 

pT range dpmjet 

QGSJET 

SIBYLL 

EPOS 

PYTHIA 

Vi = 7 TeV. 






[GeV] 

3.06 

11-04 

2.1 

LHC 

8.185 







[0.00, 0.05] 

1.24 

1.35 

0.52 

0.98 

1.13 

pT range 

DPMJET 

QGSJET 

SIBYLL 

EPOS 

PYTHIA 

[0.05, 0.10] 

1.25 

1.19 

0.49 

0.94 

1.14 

[GeV] 

3.06 

11-04 

2.1 

LHC 

8.185 

[0.10, 0.15] 

1.15 

0.91 

0.44 

0.74 

1.05 

[0.05, 0.10] 

1.96 

1.95 

1.00 

1.56 

1.92 

[0.15, 0.20] 

1.29 

0.80 

0.48 

0.83 

1.14 

[0.10, 0.15] 

1.45 

1.27 

0.69 

1.14 

1.40 

[0.20, 0.25] 

1.47 

0.74 

0.52 

0.99 

1.19 

[0.15, 0.20] 

1.27 

0.96 

0.58 

0.93 

1.21 

[0.25, 0.30] 

2.15 

0.71 

0.66 

0.80 

1.57 

[0.20, 0.25] 

1.37 

0.87 

0.62 

0.90 

1.31 

[0.30, 0.35] 

3.75 

0.57 

0.89 

0.94 

2.45 

[0.25, 0.30] 

1.39 

0.73 

0.62 

1.00 

1.28 

[0.35, 0.40] 

7.32 

0.54 

1.29 

1.38 

4.14 

[0.30, 0.35] 

1.55 

0.70 

0.69 

1.00 

1.36 

[0.40, 0.45] 

20.42 

0.64 

2.43 

3.26 

8.97 

[0.35, 0.40] 

2.20 

0.80 

0.99 

1.07 

1.86 







[0.40, 0.45] 

3.31 

0.80 

1.41 

1.21 

2.59 







[0.45, 0.50] 

4.47 

0.68 

1.78 

1.42 

3.23 







[0.50, 0.55] 

6.99 

0.67 

2.60 

1.90 

4.56 







[0.55, 0.60] 

7.84 

0.51 

2.68 

1.77 

4.50 







[0.60, 0.65] 

15.60 

0.66 

4.77 

3.27 

8.00 







[0.65, 0.70] 

14.14 

0.46 

3.70 

3.58 

6.31 







[0.70, 0.75] 

36.46 

0.79 

8.05 

14.00 

13.20 

TABLE XXIII. 

Ratio of 

7r° production rate of MC simulation 

[0.75, 0.80] 

18.03 

0.33 

2.97 

11.88 

8.13 

to data in the rapidity range 9.6 < y < 9.8 inp + p collisions 


at ^/s = 7 TeV. 

pT range dpmjet 

QGSJET 

SIBYLL 

EPOS 

PYTHIA 

[GeV] 

3.06 

11-04 

2.1 

LHC 

8.185 

[0.00, 0.05] 

1.00 

1.02 

0.35 

0.70 

0.88 

[0.05, 0.10] 

1.15 

1.03 

0.42 

0.72 

1.04 

[0.10, 0.15] 

1.28 

0.87 

0.44 

0.75 

1.10 

[0.15, 0.20] 

1.46 

0.78 

0.45 

0.97 

1.14 

[0.20, 0.25] 

2.09 

0.67 

0.53 

0.67 

1.46 

[0.25, 0.30] 

5.93 

0.72 

1.02 

1.25 

3.78 

[0.30, 0.35] 

10.39 

0.50 

1.14 

1.60 

5.78 

[0.35, 0.40] 

16.08 

0.31 

1.01 

3.17 

7.07 
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TABLE XXIV. Ratio of tt® production rate of MC simulation 


to data in the rapidity range 9.8 < y < 10.0 
at = 7 TeV. 

in p + p collisions 

pT range 
[GeV] 

DPMJET 

3.06 

QGSJET 

11-04 

SIBYLL 

2.1 

EPOS 

LHC 

PYTHIA 

8.185 

[0.00, 0.05] 

1.20 

1.13 

0.40 

0.62 

1.03 

[0.05, 0.10] 

1.38 

1.03 

0.43 

0.76 

1.14 

[0.10, 0.15] 

1.54 

0.89 

0.43 

1.03 

1.17 

[0.15, 0.20] 

2.26 

0.75 

0.49 

0.67 

1.54 

[0.20, 0.25] 

6.10 

0.65 

0.77 

1.11 

3.99 

[0.25, 0.30] 

12.79 

0.43 

0.91 

2.01 

7.39 


TABLE XXV. Ratio of tt® production rate of MC simulation 
to data in the rapidity range 10.0 < y < 10.2 in p+p collisions 
at -^/s = 7 TeV. 

pT range 
[GeV] 

DPMJET 

3.06 

QGSJET 

11-04 

SIBYLL 

2.1 

EPOS 

LHC 

PYTHIA 

8.185 

[0.00, 0.05] 

1.57 

1.14 

0.42 

1.06 

1.20 

[0.05, 0.10] 

1.73 

1.10 

0.44 

1.11 

1.24 

[0.10, 0.15] 

2.74 

0.95 

0.51 

0.77 

1.83 

[0.15, 0.20] 

9.16 

0.97 

0.90 

1.55 

6.16 

[0.20, 0.25] 

10.38 

0.31 

0.55 

1.74 

6.50 


TABLE XXVI. Ratio of tt® production rate of MC simulation 
to data in the rapidity range 10.2 < y < 10.4 in p+p collisions 
at yfs = 7 TeV. 

pT range 
[GeV] 

DPMJET 

3.06 

QGSJET 

11-04 

SIBYLL 

2.1 

EPOS 

LHC 

PYTHIA 

8.185 

[0.00, 0.05] 

2.22 

1.29 

0.45 

0.88 

1.42 

[0.05, 0.10] 

3.11 

1.05 

0.49 

0.75 

1.99 

[0.10, 0.15] 

8.51 

0.95 

0.71 

1.35 

5.94 

[0.15, 0.20] 

52.42 

1.64 

2.39 

9.00 

37.20 


TABLE XXVII. Ratio of tt® production rate of MC simulation 
to data in the rapidity range 10.4 < y < 10.6 in p+p collisions 
at = 7 TeV. 

pT range 
[GeV] 

DPMJET 

3.06 

QGSJET 

11-04 

SIBYLL 

2.1 

EPOS 

LHC 

PYTHIA 

8.185 

[0.00, 0.05] 

5.60 

1.17 

0.56 

0.97 

3.56 

[0.05, 0.10] 

9.64 

1.05 

0.64 

1.48 

7.16 

[0.10, 0.15] 

37.19 

1.39 

1.24 

7.51 

27.39 

[0.15, 0.20] 

19.23 

0.00 

0.00 

11.78 

0.00 


TABLE XXVIII. Ratio of tt® production rate of MC simula¬ 
tion to data in the rapidity range 10.6 < y < 10.8 in p + p 
collisions at ^/s — 7 TeV. 

pT range 
[GeV] 

DPMJET 

3.06 

QGSJET 

11-04 

SIBYLL 

2.1 

EPOS 

LHC 

PYTHIA 

8.185 

[0.00, 0.05] 

14.17 

0.75 

0.40 

2.06 

12.35 

[0.05, 0.10] 

17.43 

0.32 

0.50 

4.31 

11.74 

[0.10, 0.15] 

7.14 

0.00 

0.00 

2.16 

0.54 
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TABLE XXIX. Ratio of tt® production rate of MC simulation 
to data in the px range 0.0 < px < 0.2 GeV in p + p collisions 
at -s/s = 7 TeV. 


Pz range 
[GeV] 

DPMJET 

3.06 

QGSJET 

11-04 

SIBYLL 

2.1 

EPOS 

LHC 

PYTHIA 

8.185 

[ 875, 1000] 

1.20 

1.10 

0.52 

0.93 

1.12 

[1000, 1125] 

1.12 

1.00 

0.45 

0.82 

1.04 

[1125, 1250] 

0.97 

0.83 

0.37 

0.65 

0.88 

[1250, 1375] 

1.09 

0.88 

0.40 

0.61 

0.98 

[1375, 1500] 

1.20 

0.89 

0.41 

0.58 

1.03 

[1500, 1625] 

1.28 

0.88 

0.41 

0.76 

1.07 

[1625, 1750] 

1.35 

0.87 

0.41 

1.20 

1.07 

[1750, 1875] 

1.48 

0.92 

0.41 

1.04 

1.09 

[1875, 2000] 

1.69 

1.00 

0.42 

0.82 

1.17 

[2000, 2125] 

2.00 

0.99 

0.44 

0.67 

1.34 

[2125, 2250] 

2.56 

0.98 

0.51 

0.67 

1.65 

[2250, 2375] 

2.88 

0.70 

0.43 

0.63 

1.82 

[2375, 2500] 

5.17 

1.03 

0.63 

1.00 

3.27 

[2500, 2625] 

6.87 

0.98 

0.65 

1.25 

4.50 

[2625, 2750] 

8.61 

0.92 

0.64 

1.35 

6.24 

[2750, 2875] 

11.68 

0.81 

0.65 

1.48 

8.78 

[2875, 3000] 

16.87 

0.87 

0.69 

2.10 

13.53 

[3000, 3125] 

16.61 

0.45 

0.60 

2.15 

14.06 

[3125, 3250] 

21.18 

0.58 

0.66 

4.61 

19.13 

[3250, 3375] 

44.10 

0.27 

0.66 

17.76 

10.21 

[3375, 3500] 

20.03 

0.00 

0.05 

11.78 

0.99 


TABLE XXX. Ratio of tt® production rate of MC simulation 
to data in the px range 0.2 < px < 0.4 GeV in p + p collisions 
at -s/s = 7 TeV. 


Pz range 

[GeV] 

DPMJET 

3.06 

QGSJET 

11-04 

SIBYLL 

2.1 

EPOS 

LHC 

PYTHIA 

8.185 

[1000, 1125] 

1.39 

0.81 

0.72 

0.98 

1.37 

[1125, 1250] 

1.36 

0.75 

0.67 

0.93 

1.31 

[1250, 1375] 

1.39 

0.73 

0.66 

0.95 

1.30 

[1375, 1500] 

1.41 

0.70 

0.63 

0.99 

1.29 

[1500, 1625] 

1.50 

0.72 

0.65 

1.05 

1.32 

[1625, 1750] 

1.50 

0.68 

0.61 

0.95 

1.25 

[1750, 1875] 

1.64 

0.71 

0.62 

0.85 

1.31 

[1875, 2000] 

1.87 

0.71 

0.65 

0.76 

1.41 

[2000, 2125] 

2.12 

0.63 

0.65 

0.69 

1.54 

[2125, 2250] 

2.92 

0.62 

0.79 

0.80 

2.02 

[2250, 2375] 

3.89 

0.59 

0.86 

0.95 

2.57 

[2375, 2500] 

4.93 

0.58 

0.93 

1.07 

3.11 

[2500, 2625] 

10.24 

0.92 

1.59 

2.04 

6.40 

[2625, 2750] 

12.52 

0.81 

1.57 

1.92 

7.62 

[2750, 2875] 

18.24 

0.73 

1.75 

2.33 

10.62 

[2875, 3000] 

22.25 

0.68 

1.61 

2.99 

13.15 

[3000, 3125] 

30.98 

0.55 

1.45 

4.81 

16.99 

[3125, 3250] 

36.45 

0.46 

1.81 

9.28 

21.07 

[3250, 3375] 

39.91 

0.25 

1.05 

13.29 

6.50 

[3375, 3500] 

39.63 

0.05 

0.00 

15.28 

2.05 


TABLE XXXI. Ratio of tt® production rate of MC simulation 
to data in the px range 0.4 < px < 0.6 GeV in p + p collisions 
at = 7 TeV. 


Pz range 

[GeV] 

DPMJET 

3.06 

QGSJET 

11-04 

SIBYLL 

2.1 

EPOS 

LHC 

PYTHIA 

8.185 

[1500, 1625] 

2.66 

0.89 

1.65 

1.35 

2.46 

[1625, 1750] 

2.41 

0.72 

1.41 

1.11 

2.14 

[1750, 1875] 

2.27 

0.60 

1.25 

0.96 

1.94 

[1875, 2000] 

2.82 

0.63 

1.44 

1.10 

2.30 

[2000, 2125] 

3.51 

0.64 

1.64 

1.21 

2.75 

[2125, 2250] 

3.42 

0.48 

1.45 

1.10 

2.51 

[2250, 2375] 

4.90 

0.56 

1.85 

1.45 

3.36 

[2375, 2500] 

5.95 

0.53 

1.92 

1.52 

3.74 

[2500, 2625] 

7.72 

0.55 

2.16 

1.61 

4.48 

[2625, 2750] 

13.14 

0.70 

3.03 

2.22 

6.82 

[2750, 2875] 

20.33 

0.73 

3.66 

3.38 

9.43 

[2875, 3000] 

19.78 

0.52 

2.69 

3.09 

8.08 



37 


TABLE XXXII. Ratio of tt® production rate of MC simulation 
to data in the px range 0.6 < px < 0.8 GeV in p + p collisions 
at -s/s = 7 TeV. 


Pz range 
[GeV] 

DPMJET 

3.06 

QGSJET 

11-04 

SIBYLL 

2.1 

EPOS 

LHC 

PYTHIA 

8.185 

[2000, 2125] 

3.18 

0.47 

2.32 

1.25 

2.80 

[2125, 2250] 

4.30 

0.53 

2.85 

1.52 

3.60 

[2250, 2375] 

4.50 

0.46 

2.81 

1.44 

3.51 

[2375, 2500] 

6.52 

0.55 

3.71 

1.83 

4.98 

[2500, 2625] 

7.37 

0.50 

3.69 

1.93 

5.06 

[2625, 2750] 

8.01 

0.45 

3.51 

1.85 

5.05 

[2750, 2875] 

9.74 

0.37 

3.53 

2.03 

5.40 

[2875, 3000] 

15.66 

0.46 

4.41 

3.52 

7.32 

[3000, 3125] 

22.78 

0.50 

4.98 

5.70 

8.80 

[3125, 3250] 

18.98 

0.39 

3.92 

7.66 

7.07 

[3250, 3375] 

10.02 

0.13 

0.89 

5.85 

2.73 


TABLE XXXIII. Ratio of tt® production rate of MC simula¬ 
tion to data in the px range 0.8 < px < 1.0 GeV in p + p 
collisions at ■•/s = 7 TeV. 


Pz range 

[GeV] 

DPMJET 

3.06 

QGSJET 

11-04 

SIBYLL 

2.1 

EPOS 

LHC 

PYTHIA 

8.185 

[2250, 2375] 

4.75 

0.45 

4.34 

1.72 

4.54 

[2375, 2500] 

6.50 

0.55 

5.55 

2.45 

5.91 

[2500, 2625] 

8.07 

0.57 

6.18 

2.67 

6.95 

[2625, 2750] 

11.27 

0.59 

7.67 

3.30 

9.16 

[2750, 2875] 

9.22 

0.35 

5.65 

2.56 

7.02 

[2875, 3000] 

22.52 

0.65 

10.80 

7.06 

15.46 

[3000, 3125] 

16.96 

0.32 

6.48 

6.12 

10.31 

[3125, 3250] 

46.87 

0.90 

15.67 

23.77 

26.10 

[3250, 3375] 

33.52 

0.43 

7.16 

30.40 

17.58 
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TABLE XXXIV. Ratio of tt® production rate of MC simu¬ 
lation to data in the rapidity range 8.8 < y < 9.0 in p + p 
collisions at -y/s = 2.76 TeV. 


pT range 
[GeV] 

DPMJET 

3.06 

QGSJET 

11-04 

SIBYLL 

2.1 

EPOS 

LHC 

PYTHIA 

8.185 

[0.10, 0.15] 

1.67 

1.13 

0.61 

1.24 

1.45 

[0.15, 0.20] 

2.04 

0.93 

0.62 

0.87 

1.59 

[0.20, 0.25] 

6.68 

1.23 

1.36 

1.57 

4.90 

[0.25, 0.30] 

20.32 

1.32 

2.43 

3.07 

14.56 


TABLE XXXVIII. Ratio of production rate of MC simu¬ 
lation to data in the rapidity range 9.6 < y < 9.8 in p -\- p 
collisions at y/s = 2.76 TeV. 


pT range 
[GeV] 

DPMJET 

3.06 

QGSJET 

11-04 

SIBYLL 

2.1 

EPOS 

LHC 

PYTHIA 

8.185 

[0.00, 0.05] 

4.63 

0.51 

0.38 

0.56 

4.35 

[0.05, 0.10] 

8.44 

0.48 

0.36 

1.06 

7.80 

[0.10, 0.15] 

5.46 

0.06 

0.13 

1.66 

1.81 


TABLE XXXV. Ratio of tt® production rate of MC simulation 
to data in the rapidity range 9.0 < y < 9.2 in p -\-p collisions 
at y/s = 2.76 TeV. 


pT range 
[GeV] 

DPMJET 

3.06 

QGSJET 

11-04 

SIBYLL 

2.1 

EPOS 

LHC 

PYTHIA 

8.185 

[0.05, 0.10] 

1.47 

1.08 

0.49 

1.07 

1.20 

[0.10, 0.15] 

1.85 

0.91 

0.51 

0.77 

1.40 

[0.15, 0.20] 

5.73 

1.01 

0.93 

1.18 

4.16 

[0.20, 0.25] 

34.50 

2.00 

3.01 

4.31 

26.06 


TABLE XXXVI. Ratio of tt® production rate of MC simu¬ 
lation to data in the rapidity range 9.2 < y < 9.4 in p + p 
collisions at -/s = 2.76 TeV. 


pT range 
[GeV] 

DPMJET 

3.06 

QGSJET 

11-04 

SIBYLL 

2.1 

EPOS 

LHC 

PYTHIA 

8.185 

[0.00, 0.05] 

1.22 

0.89 

0.34 

0.72 

0.91 

[0.05, 0.10] 

2.00 

0.98 

0.46 

0.69 

1.41 

[0.10, 0.15] 

3.52 

0.64 

0.46 

0.65 

2.54 

[0.15, 0.20] 

10.45 

0.63 

0.73 

1.43 

8.89 

[0.20, 0.25] 

15.42 

0.34 

0.33 

4.37 

7.87 


TABLE XXXVII. Ratio of tt® production rate of MC simu¬ 
lation to data in the rapidity range 9.4 < y < 9.6 in p + p 
collisions at y/s = 2.76 TeV. 


pT range 
[GeV] 

DPMJET 

3.06 

QGSJET 

11-04 

SIBYLL 

2.1 

EPOS 

LHC 

PYTHIA 

8.185 

[0.00, 0.05] 

2.12 

0.72 

0.31 

0.47 

1.36 

[0.05, 0.10] 

3.65 

0.70 

0.39 

0.59 

2.86 

[0.10, 0.15] 

11.29 

0.68 

0.69 

1.48 

9.99 

[0.15, 0.20] 

9.89 

0.19 

0.18 

3.04 

4.37 
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TABLE XXXIX. Ratio of tt® production rate of MC simula¬ 
tion to data in the pr range 0.0 < pr < 0.2 GeV in p -|- p 
collisions at = 2.76 TeV. 


Pz range 
[GeV] 

DPMJET 

3.06 

QGSJET 

11-04 

SIBYLL 

2.1 

EPOS 

LHC 

PYTHIA 

8.185 

[ 625, 750] 

1.34 

0.94 

0.46 

1.05 

1.11 

[ 750, 875] 

1.66 

0.89 

0.45 

0.63 

1.22 

[ 875, 1000] 

2.63 

0.65 

0.42 

0.56 

1.79 

[1000, 1125] 

5.37 

0.70 

0.53 

0.83 

4.42 

[1125, 1250] 

12.08 

0.60 

0.63 

1.30 

11.26 

[1250, 1380] 

18.30 

0.41 

0.43 

4.88 

10.87 
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TABLE XL. Ratio of production rate of MC simulation to 
data in the px range 0.2 < pt < 0.4 GeV in p + p collisions at 
Vi = 2.76TeV. 


Pz range 
[GeV] 

DPMJET 

3.06 

QGSJET 

11-04 

SIBYLL 

2.1 

EPOS 

LHC 

PYTHIA 

8.185 

[ 750, 875] 

2.56 

0.84 

0.93 

0.93 

2.03 

[ 875, 1000] 

5.10 

0.79 

1.24 

1.28 

3.66 

[1000, 1125] 

11.69 

0.87 

1.73 

1.82 

7.99 

[1125, 1250] 

138.29 

4.19 

10.41 

16.28 

90.57 

[1250, 1380] 

44.14 

0.45 

1.61 

13.18 

15.99 
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TABLE XLI. Ratio of tt® production rate of MC simulation 
to data in the rapidity range —8.8 > piab > —9.0 in p + Pb 
collisions at y^SNN = 5.02 TeV. 

pT range 
[GeV] 

DPMJET 

3.06 

QGSJET 

11-04 

EPOS 

LHC 

[0.10, 

0.15] 

1.39 

1.41 

1.15 

[0.15, 

0.20] 

0.99 

0.92 

0.85 

[0.20, 

0.25] 

0.90 

0.77 

0.77 

[0.25, 

0.30] 

0.83 

0.62 

0.66 

[0.30, 

0.35] 

0.94 

0.60 

0.70 

[0.35, 

0.40] 

0.91 

0.53 

0.63 

[0.40, 

0.45] 

1.03 

0.55 

0.69 

[0.45, 

0.50] 

1.35 

0.65 

0.85 

[0.50, 

0.55] 

1.28 

0.56 

0.79 

[0.55, 

0.60] 

1.62 

0.58 

0.93 

[0.60, 

0.65] 

1.89 

0.60 

1.07 

[0.65, 

0.70] 

2.11 

0.58 

1.15 

[0.70, 

0.75] 

2.03 

0.51 

1.11 

[0.75, 

0.80] 

3.44 

0.79 

1.77 

[0.80, 

0.85] 

5.00 

1.05 

2.68 

[0.85, 

0.90] 

2.83 

0.56 

1.75 

[0.90, 

0.95] 

10.56 

1.88 

6.44 

[0.95, 

1.00] 

29.43 

4.41 

23.12 


TABLE XLII. Ratio of rr® production rate of MC simulation 
to data in the rapidity range —9.0 > piab > —9.2 in p + Pb 
collisions at y^SNN = 5.02 TeV. 

pT range 

DPMJET 

QGSJET 

EPOS 

[GeV] 

3.06 

11-04 

LHC 

[0.05, 0.10] 

1.41 

1.57 

1.15 

[0.10, 0.15] 

1.05 

1.03 

0.87 

[0.15, 0.20] 

0.86 

0.78 

0.73 

[0.20, 0.25] 

0.85 

0.71 

0.71 

[0.25, 0.30] 

0.98 

0.68 

0.73 

[0.30, 0.35] 

0.96 

0.57 

0.65 

[0.35, 0.40] 

1.14 

0.64 

0.74 

[0.40, 0.45] 

1.30 

0.65 

0.79 

[0.45, 0.50] 

1.43 

0.62 

0.81 

[0.50, 0.55] 

1.51 

0.49 

0.77 

[0.55, 0.60] 

2.80 

0.75 

1.32 

[0.60, 0.65] 

3.33 

0.73 

1.52 

[0.65, 0.70] 

2.59 

0.50 

1.13 

[0.70, 0.75] 

6.67 

1.14 

2.88 

[0.75, 0.80] 

4.88 

0.87 

2.48 


TABLE XLIII. Ratio of tt® production rate of MC simulation 
to data in the rapidity range —9.2 > j/iab > —9.4 in p + Pb 
collisions at Asnn = 5.02 TeV. 


pT range 
[GeV] 

DPMJET 

3.06 

QGSJET 

11-04 

EPOS 

LHC 

[0.00, 0.05] 

1.39 

1.68 

1.12 

[0.05, 0.10] 

1.08 

1.16 

0.88 

[0.10, 0.15] 

0.92 

0.87 

0.75 

[0.15, 0.20] 

0.79 

0.70 

0.67 

[0.20, 0.25] 

0.88 

0.69 

0.68 

[0.25, 0.30] 

1.09 

0.69 

0.75 

[0.30, 0.35] 

1.16 

0.67 

0.75 

[0.35, 0.40] 

1.32 

0.66 

0.75 

[0.40, 0.45] 

1.41 

0.56 

0.71 

[0.45, 0.50] 

1.80 

0.49 

0.76 

[0.50, 0.55] 

5.02 

1.09 

1.86 

[0.55, 0.60] 

4.57 

0.81 

1.63 


TABLE XLIV. Ratio of tt® production rate of MC simulation 
to data in the rapidity range —9.4 > j/iab > —9.6 in p + Pb 
collisions at y^SNN = 5.02 TeV. 

pT range 

DPMJET 

QGSJET 

EPOS 

[GeV] 

3.06 

11-04 

LHC 

[0.00, 0.05] 

1.12 

1.27 

0.90 

[0.05, 0.10] 

0.98 

0.99 

0.78 

[0.10, 0.15] 

0.86 

0.79 

0.71 

[0.15, 0.20] 

0.85 

0.71 

0.67 

[0.20, 0.25] 

1.02 

0.67 

0.70 

[0.25, 0.30] 

1.14 

0.68 

0.71 

[0.30, 0.35] 

0.95 

0.46 

0.50 

[0.35, 0.40] 

1.27 

0.42 

0.54 

[0.40, 0.45] 

2.64 

0.59 

0.89 


TABLE XLV. Ratio of tt® production rate of MC simulation 
to data in the rapidity range —9.6 > j/iab > —9.8 in p + Pb 
collisions at y^SNN = 5.02 TeV. 

pT range 

DPMJET 

QGSJET 

EPOS 

[GeV] 

3.06 

11-04 

LHC 

[0.00, 0.05] 

0.96 

1.02 

0.78 

[0.05, 0.10] 

0.88 

0.84 

0.71 

[0.10, 0.15] 

0.92 

0.80 

0.73 

[0.15, 0.20] 

1.04 

0.73 

0.72 

[0.20, 0.25] 

1.19 

0.73 

0.72 

[0.25, 0.30] 

1.36 

0.62 

0.64 

[0.30, 0.35] 

2.73 

0.77 

0.96 

[0.35, 0.40] 

9.27 

1.71 

2.57 
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TABLE XLVI 
to data in the 
collisions at 

. Ratio of 7r° production rate of MC simulation 
rapidity range —9.8 > j/iab > —10.0 in p + Pb 
snn = 5.02 TeV. 

pT range 

DPMJET 

QGSJET 

EPOS 

[GeV] 

3.06 

11-04 

LHC 

[0.00, 0.05] 

0.86 

0.85 

0.68 

[0.05, 0.10] 

0.92 

0.83 

0.71 

[0.10, 0.15] 

1.18 

0.89 

0.83 

[0.15, 0.20] 

1.27 

0.80 

0.75 

[0.20, 0.25] 

1.23 

0.54 

0.53 

[0.25, 0.30] 

2.98 

0.73 

0.90 


TABLE XLVII. Ratio of tt® production rate of MC simulation 
to data in the rapidity range —10.0 > piab > —10.2 in p + Pb 
collisions at y^SNN = 5.02 TeV. 

pT range 

DPMJET 

QGSJET 

EPOS 

[GeV] 

3.06 

11-04 

LHC 

[0.00, 0.05] 

1.00 

0.90 

0.74 

[0.05, 0.10] 

1.14 

0.91 

0.81 

[0.10, 0.15] 

1.77 

1.15 

1.00 

[0.15, 0.20] 

1.28 

0.55 

0.52 

[0.20, 0.25] 

4.66 

1.06 

1.31 


TABLE XLVIII. Ratio of production rate of MC simulation 
to data in the rapidity range —10.2 > j/iab > —10.4 in p + Pb 
collisions at y^SNN = 5.02 TeV. 

pT range 

DPMJET 

QGSJET 

EPOS 

[GeV] 

3.06 

11-04 

LHC 

[0.00, 0.05] 

1.73 

1.43 

1.25 

[0.05, 0.10] 

1.52 

1.06 

0.85 

[0.10, 0.15] 

1.56 

0.68 

0.60 

[0.15, 0.20] 

4.02 

0.89 

1.04 


TABLE XLIX. Ratio of production rate of MC simulation 
to data in the rapidity range —10.4 > piab > —10.6 in p + Pb 
collisions at = 5.02 TeV. 

pT range 

DPMJET 

QGSJET 

EPOS 

[GeV] 

3.06 

11-04 

LHC 

[0.00, 0.05] 

1.40 

0.92 

0.75 

[0.05, 0.10] 

2.68 

1.18 

0.93 

[0.10, 0.15] 

5.58 

1.32 

1.45 

[0.15, 0.20] 

1.30 

0.15 

0.39 


TABLE L. Ratio of tt® production rate of MC simulation to 
data in the rapidity range —10.6 > piab > —10.8 in p + Pb 
collisions at y^SNN = 5.02 TeV. 

pT range 

DPMJET 

QGSJET 

EPOS 

[GeV] 

3.06 

11-04 

LHC 

[0.00, 0.05] 

2.17 

0.97 

0.98 

[0.05, 0.10] 

5.77 

1.52 

1.61 

[0.10, 0.15] 

4.60 

0.65 

1.47 
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TABLE LI. Ratio of 7r° production rate of MC simulation to 
data in the Pt range 0.0 < pr < 0.2 GeV in p + Pb collisions 
at Vsnn = 5.02 TeV. 


Pz range 

[GeV] 

DPMJET 

3.06 

QGSJET 

11-04 

EPOS 

LHC 

[ 875, fOOO] 

0.97 

0.97 

0.80 

[1000, 1125] 

0.88 

0.86 

0.73 

[1125, 1250] 

0.78 

0.74 

0.65 

[1250, 1375] 

0.80 

0.74 

0.65 

[1375, 1500] 

0.85 

0.77 

0.69 

[1500, 1625] 

0.93 

0.81 

0.70 

[1625, 1750] 

1.11 

0.92 

0.79 

[1750, 1875] 

1.20 

0.92 

0.84 

[1875, 2000] 

1.55 

1.13 

1.14 

[2000, 2125] 

1.27 

0.93 

0.86 

[2125, 2250] 

1.77 

1.26 

1.04 

[2250, 2375] 

1.61 

1.09 

0.86 

[2375, 2500] 

1.38 

0.83 

0.66 

[2500, 2625] 

1.48 

0.78 

0.66 

[2625, 2750] 

1.77 

0.79 

0.73 

[2750, 2875] 

1.73 

0.69 

0.66 

[2875, 3000] 

1.77 

0.63 

0.60 

[3000, 3125] 

2.23 

0.71 

0.70 

[3125, 3250] 

5.20 

1.47 

1.54 

[3250, 3375] 

5.92 

1.40 

1.66 

[3375, 3500] 

6.44 

1.57 

1.65 

[3500, 3625] 

10.06 

1.75 

2.85 

[3625, 3750] 

8.63 

1.53 

3.11 

[3750, 3875] 

25.57 

1.70 

10.48 

[3875, 4000] 

10.74 

0.70 

7.90 


TABLE LII. Ratio of production rate of MC simulation to 
data in the pr range 0.2 < px < 0.4 GeV in p + Pb collisions 
at y^SNN = 5.02 TeV. 


Pz range 

[GeV] 

DPMJET 

3.06 

QGSJET 

11-04 

EPOS 

LHC 

[ 875, 1000] 

0.88 

0.68 

0.73 

[1000, 1125] 

0.76 

0.58 

0.62 

[1125, 1250] 

0.77 

0.57 

0.61 

[1250, 1375] 

0.79 

0.56 

0.60 

[1375, 1500] 

0.91 

0.61 

0.65 

[1500, 1625] 

1.02 

0.65 

0.71 

[1625, 1750] 

0.94 

0.58 

0.64 

[1750, 1875] 

1.02 

0.61 

0.68 

[1875, 2000] 

1.23 

0.73 

0.81 

[2000, 2125] 

1.30 

0.76 

0.78 

[2125, 2250] 

1.28 

0.72 

0.72 

[2250, 2375] 

1.31 

0.68 

0.67 

[2375, 2500] 

1.24 

0.57 

0.59 

[2500, 2625] 

1.16 

0.46 

0.52 

[2625, 2750] 

1.89 

0.61 

0.72 

[2750, 2875] 

2.25 

0.63 

0.78 

[2875, 3000] 

2.07 

0.53 

0.66 

[3000, 3125] 

3.27 

0.68 

0.88 

[3125, 3250] 

5.21 

0.92 

1.25 

[3250, 3375] 

4.39 

0.68 

1.03 

[3375, 3500] 

7.73 

0.99 

1.79 

[3500, 3625] 

4.46 

0.36 

1.02 

[3625, 3750] 

2.53 

0.23 

0.93 

[3750, 3875] 

1.81 

0.11 

0.58 

[3875, 4000] 

1.07 

0.05 

0.74 


TABLE LIII. Ratio of tt® production rate of MC simulation to 
data in the px range 0.4 < px < 0.6 GeV in p + Pb collisions 
at snn = 5.02 TeV. 


Pz range 
[GeV] 

DPMJET 

3.06 

QGSJET 

11-04 

EPOS 

LHC 

[1375, 1500] 

0.89 

0.45 

0.62 

[1500, 1625] 

1.23 

0.60 

0.83 

[1625, 1750] 

1.05 

0.51 

0.69 

[1750, 1875] 

1.11 

0.52 

0.70 

[1875, 2000] 

1.24 

0.57 

0.77 

[2000, 2125] 

1.30 

0.58 

0.77 

[2125, 2250] 

1.36 

0.58 

0.78 

[2250, 2375] 

1.39 

0.53 

0.73 

[2375, 2500] 

2.19 

0.75 

1.07 

[2500, 2625] 

1.64 

0.47 

0.76 

[2625, 2750] 

2.28 

0.61 

0.97 

[2750, 2875] 

3.45 

0.81 

1.32 

[2875, 3000] 

1.72 

0.36 

0.61 
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TABLE LIV. Ratio of tt® production rate of MC simulation to 
data in the Pt range 0.6 < pr < 0.8 GeV in p + Pb collisions 
at Vsnn = 5.02 TeV. 


Pz range 

[GeV] 

DPMJET 

3.06 

QGSJET 

11-04 

EPOS 

LHC 

[1750, 1875] 

1.35 

0.49 

0.93 

[1875, 2000] 

1.25 

0.44 

0.84 

[2000, 2125] 

1.08 

0.37 

0.71 

[2125, 2250] 

1.44 

0.46 

0.89 

[2250, 2375] 

1.42 

0.44 

0.86 

[2375, 2500] 

1.52 

0.43 

0.86 

[2500, 2625] 

1.86 

0.48 

1.02 

[2625, 2750] 

2.79 

0.68 

1.43 

[2750, 2875] 

2.38 

0.52 

1.14 

[2875, 3000] 

4.95 

1.05 

2.26 

[3000, 3125] 

8.90 

1.73 

3.49 

[3125, 3250] 

2.12 

0.38 

0.83 

[3250, 3375] 

3.24 

0.52 

1.29 

[3375, 3500] 

5.59 

0.85 

2.10 


TABLE LV. Ratio of tt® production rate of MC simulation to 
data in the pr range 0.8 < px < 1.0 GeV in p + Pb collisions 
at .Vsnn = 5.02 TeV. 


Pz range 

[GeV] 

DPMJET 

3.06 

QGSJET 

11-04 

EPOS 

LHC 

[2375, 2500] 

1.68 

0.47 

1.20 

[2500, 2625] 

1.64 

0.41 

1.09 

[2625, 2750] 

2.18 

0.50 

1.44 

[2750, 2875] 

1.86 

0.42 

1.19 

[2875, 3000] 

2.14 

0.46 

1.24 

[3000, 3125] 

3.36 

0.69 

1.93 

[3125, 3250] 

9.14 

1.78 

5.46 

[3250, 3375] 

15.02 

2.75 

9.19 

[3375, 3500] 

8.84 

1.35 

5.63 
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TABLE LVI. Nuclear modification factor for the 7r°s in the rapidity ranges —8.8 > j/iab > —9.0, —9.0 > j/iab > —9.2, and 
—9.2 > j/iab > —9.4 in p + Pb collisions at ^snn = 5.02 TeV. The uncertainties include the both statistical and systematic 
uncertainties. 


Pt range [GeV] 


Nuclear modification factors 



-8.8 > piab > -9.0 

-9.0 > piab > -9.2 

—9.2 > piab > —9.4 


[0.00,0.05] 

[0.05,0.10] 

[0.10,0.15] 

[0.15,0.20] 

[0.20,0.25] 

[0.25,0.30] 

[0.30,0.35] 

[0.35,0.40] 

[0.40,0.45] 

[0.45,0.50] 

[0.50,0.55] 

[0.55,0.60] 

[0.60,0.65] 

[0.65,0.70] 

[0.70,0.75] 

[0.75,0.80] 

[0.80,0.85] 

[0.85,0.90] 

[0.90,0.95] 


"-4.70x10“ 


l.llxl0“^-5.42x10“ 


1.49x10“^ 

2.11x10“^ 

2.29x10"^ 

2.28x10“^ 

1.84x10“^ 

2.03x10“^ 

1.95x10“^ 

1.59x10“^ 

2.10x10“^ 

2.06x10“^ 

2.25x10"^ 

2.74x10“^ 

4.18x10“^ 

2.69x10“^ 

1.63x10"^ 

8.59x10“^ 

1.10x10“® 


-7.67x10“ 

-9.22x10“ 

-9.86x10“ 

-7.71x10“ 

-6.09x10“ 

-7.26x10“ 

-8.44x10“ 

-8.32x10“ 

-1.24x10“ 

-1.40x10“ 

-1.73x10“ 

-2.32x10“ 

-3.84x10“ 

-2.69x10“ 

-1.63x10“ 

-8.59x10“ 

-1.10x10“ 


"+4.70x10 

^+5.42x10 

^+7.67x10 

^+9.22x10 

^+9.86x10 

^+7.71x10 

^+6.09x10 

^+7.26x10 

^+8.44x10 

^+8.32x10 

^+1.24x10 

^+1.40x10 

^+1.73x10 

^+2.32x10 

^+3.84x10 

^+2.94x10 

^+2.66x10 

^+9.92x10 

"+1.16x10 


-1 
1 

-1 
1 

-1 
1 

-1 
1 

-1 
1 

-2 
2 

2.29x10“^ 

3.31x10“^ 

9.71x10“^ 


1.12x10“ 

1.42x10“ 

1.76x10“ 

1.69x10“ 

1.21x10“ 

1.32x10“ 

1.17x10“ 

1.17x10“ 

1.28x10“ 

1.67x10“ 

8.96x10“ 

9.99x10“ 


-4.92x10“ 

-5.86x10“ 

-7.39x10“ 

-7.19x10“ 

-4.44x10“ 

-4.75x10“ 

-5.32x10“ 

-6.55x10“ 

-8.46x10“ 

-1.24x10“ 

-8.25x10“ 

-9.99x10“ 

-2.29x10“ 

-3.31x10“ 

-9.71x10“ 


"+4.92x10 

"+5.86x10 

"+7.39x10“ 

"+7.19x10“ 

"+4.44x10“ 

"+4.75x10“ 

"+5.32x10“ 

"+6.55x10“ 

"+8.46x10“ 

1+1.24x10“ 

"+8.25x10“ 

"+1.04x10“ 

1+2.31x10“ 

"+1.09x10“ 

"+1.96x10“ 


-2 

“" 1.25xl0"i- 
“" 1.60xl0"i- 
“" 1.15xl0"i- 


■" 8.53x10""- 
■" 8.89x10""- 
■" 9.24x10""- 
■" 1.10x10"!- 
■" 1.12x10"!- 
■1 2.10x10""- 

■" 4.80x10""- 
-1 

-1 

-1 


5.20x10 

6.79x10 

4.61x10 

3.26x10 

4.08x10 

5.55x10 

8.05x10 

1.03x10 

2.10x10 

4.80x10 


“"+5.20x10"" 

“"+6.79x10"" 

“"+4.61x10"" 

“"+3.26x10"" 

“"+4.08x10"" 

“"+5.55x10"" 

“"+8.05x10"" 

“1+1.03x10"! 

“"+5.78x10"" 

“"+9.11x10"" 


TABLE LVII. Nuclear modification factor for the tt^s in the rapidity ranges —9.4 > j/iab > —9.6, —9.6 > piab > —9.8, and 
—9.8 > i/iab > —10.0 in p + Pb collisions at y^SNN = 5.02 TeV. The uncertainties include the both statistical and systematic 
uncertainties. 


Pt range [GeV] 


Nuclear modification factors 




-9.4 > piab > -9.6 

-9.6 > piab > -9.8 

-9.8 > piab > -10.0 



[0.00,0.05] 8.99xl0""-4.55xl0""+4.55x10"" 9.56x 10""-5.10x 10""+5.10x 10"" 1.13x 10"i-6.96x 10""+6.96x 10"" 

[0.05,0.10] 1.02X 10"i-4.60x 10""+4.60x 10"" 1.09x 10"i-5.36x 10""+5.36x 10"" 9.91x 10""-5.54x 10""+5.54x 10"" 

[0.10,0.15] 1.18xl0"i-5.41xl0""+5.41xl0"" 9.78x 10""-4.90x 10""+4.90x 10"" 6.62xl0""-3.11xl0""+3.11xl0"" 

[0.15,0.20] 1.15X 10"i-4.96x 10""+4.96x 10"" 7.98x 10""-3.08x 10""+3.08x 10"" 7.79x 10""-3.47x 10""+3.47x 10"" 

[0.20,0.25] 8.38X 10""-3.11x 10""+3.11x 10"" 7.81x 10""-3.51x 10""+3.51x 10"" 1.18x 10"i-7.62x 10""+7.62x 10"" 

[0.25,0.30] 8.27X 10""-3.76x 10""+3.76x 10"" 8.54x 10""-5.53x 10""+5.53x 10"" 4.45x 10""-4.45x 10""+5.07x 10"" 

[0.30,0.35] 1.37X 10"i-7.89x 10""+7.89x 10"" 4.00x 10""-4.00x 10""+4.21x 10"" 

[0.35,0.40] 1.24xl0"i-9.88xl0""+9.88x10"" l.lOx 10"®-1.10x 10"®+6.11x 10"" 

[0.40,0.45] 6.08X 10""-6.08x 10""+7.41x 10"" 
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TABLE LVIII. Nuclear modification factor for the 7r°s in the rapidity ranges —10.0 > t/iab > —10.2, —10.2 > yiab > —10.4, and 
— 10.4 > 2/iab > —10.6 in p + Pb collisions at ^'snn = 5.02 TeV. The uncertainties include the both statistical and systematic 
uncertainties. 


Pt range [GeV] 

-10.0 > piab > -10.2 

[0.00,0.05] 7.83X 10~''-5.73x 10'^+5.73x lO’"' 

[0.05,0.10] 7.16X 10“^-3.70x 10~2+3.70 x 10"^ 

[0.10,0.15] 4.81 X 10“^-2.71x 10"^+2.71 x 10"^ 

[0.15,0.20] 1.42X 10“^-9.45x 10”^+9.45x 10"^ 

[0.20,0.25] 1.92X 10“^-1.92x 10"^+4.11 X10"^ 


Nuclear modification factors 
-10.2 > piab > -10.4 
1.89xl0-''-1.89xl0“''+2.61xl0~^ 
6.42X 10"^-2.98x 10“^+2.98x 10"^ 
1.38X 10~^-9.05x 10“^+9.05x 10"^ 
4.75X 10"^-4.75x 10“^+5.26x 10"^ 


-10.4 > piab > -10.6 
3.92X 10"^-3.55x 10-''+3.55x lO""" 
4.78X 10"^-4.78x 10"^+5.02x 10"^ 
1.32xl0"^-1.32xl0"^+5.04x10"^ 
2.06xl0"^-2.06xl0"^+2.38x10"^ 


TABLE LIX. Nuclear modification factor for the tt'^s in the rapidity ranges —10.6 > piab > —10.8 in p + Pb collisions at 
^snn = 5.02 TeV. The uncertainties include the both statistical and systematic uncertainties. 


Pt range ]GeV] 

Nuclear modification factors 

-10.6 > Plab > -10.8 

[0.00,0.05] 

1.44xl0"''-1.44xl0“''+1.96xl0'-^ 

[0.05,0.10] 

1.10xl0"®-1.10xl0“®+5.60x10"^ 

[0.10,0.15] 

3.00X 10"^-3.00x 10“^+4.24x 10"^ 

[0.15,0.20] 

4.75X 10~^-4.75x 10“^+5.26x 10~^ 



